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W. Pogorzelski: 


Badania teoretyczne ilości ciepła, otrzymy- 
wanych na kuli ziemskiej, z uwzględnieniem 
strat promieniowania w amosferze. 


(Z Pracowni Meteorologicznej przy Tow. Nauk. Warsz.). 


Przedstawił Wł. Gorczyński. 


l. Rozważania ogólne. 


Badania niniejsze były tematem licznych rozpraw, z któ- 
rych ważniejsze, wraz z dziełami, odnoszącemi się do przedmio- 
tu insolacyi, podane są poniżej w wykazie bibliograficznym. 

Obliczenia ilości ciepła, których będziemy dokonywali je- 
dynie przy pomocy teoretycznych rozumowań, przeprowadzimy 
najpierw w założeniu, iż całkowita ilość energii promienistej do- 
chodzi do powierzchni ziemi, lub, co na jedno wychodzi, dla jed- 
nostki powierzchni poziomej, położonej na granicy atmosfery. 
Biorąc następnie za podstawę prawo Lamberta, które jest 
jednem z prostszych wyrażeń dla strat promieniowania w atmo- 
sferze, dojdziemy do nowych wyrażeń ilości ciepła. 

Następnym przedmiotem rozważań będzie przeprowadzenie 
podobnych rachunków, jak również i krytyka poprzednich rezul- 
tatów, oparte na wynikach badań w dziedzinie spektrobolo- 
metryi. 


Jeśli przez h oznaczymy wysokość słońca, przez r odległość 
środka słońca od emt, to ilość ciepła dq, jaką otrzyma jednost- 
ka powierzchni ziemi w czasie dź, wyrazi się wzorem: 


0) dq SĘ Sinh. dt 


FER | 


Stata C będzie wyrażała ilość ciepła, otrzymaną w jednost- 
ce odległości przez jednostkę powierzchni, prostopadłej do pro- 


mieni słońca. Wielkość będzie równa ilości ciepła, otrzyma- 


nego na granicy atmosfery; jest to tak zwana stała słoneczna, 
która ulega zresztą drobnym wahaniom z powodu zmienności 7. 

Dla odległości średniej od słońca, według pomiarów W. Gor- 
czyńskiego !), wartość stałej słonecznej wynosi 2,05 kal. gr. 
cm.*/min.; w rachunkach następnych przyjmować będziemy dla sta- 
łej słonecznej wartości 2 kal. gr. 

Ilość ciepła, którą otrzyma w ciągu dnia jednostka po- 
wierzchni ziemi o danem położeniu geograficznem, obliczyć moż- 
na, całkując (1) w granicach, odpowiadających wschodowi i za- 
chodowi słońca. Ilość dzienna ciepła będzie zatem pewną funk- 
cyą długości słońca; należy jednak zauważyć, iż długość będzie 
się tu zmieniała z jednego dnia na drugi w sposób nieciągły, 
gdyż z chwilą nadejścia nocy, gdy ciepło już nie dochodzi, dłu- 
gość słońca zmienia się w dalszym ciągu. Weźmy teraz pod 
uwagę ilość ciepła, otrzymaną w dłuższym przeciągu czasu, od- 
powiadającego długościom słońca (©,, ©); jedyną ścisłą i bez- 
pośrednią metodą będzie zsumowanie ilości ciepła dla wszyst- 
kich poszczególnych dni, zawartych w okresie (©,, ©.). Oznacz- 
my więc przez f (7) natężenie promieniowania, gdzie + jest ką- 
tem godzinnym słońca, otrzymamy dla ilości ciepła wyrażenie 
następujące : 


*0 


©: WAT V 
(2) Q E f(x) d 
8-8 dl MA: 


to 


* będzie połową łuku dziennego słońca. 

Powyższy rachunek upraszcza się przez wprowadzenie po- 
jęcia ilości średniej ciepła, sprowadza się bowiem do całkowa- 
nia. Metoda ta przedstawioną będzie w rozdziale następnym. 


1) W. Gorczyński: O wartości stałej słonecznej, według pomiarów 
dokonanych sposobem ekranów absorbcyjnych w Ursynowie. (Sprawozdania 
Tow. Nauk. Warsz., 1909 r.). 


OE 


2. llości ciepła dla równika i bieguna. 


Jak wynika z poprzedniego, ilość dzienna ciepła będzie: 


(1) os = fSinn.at 
Oznaczywszy przez 4 — szerokość geograficzną danego 
miejsca, przez à — zboczenie słońca, mamy z trójkąta parala- 


ktycznego: 

(2) Sinh = Sing. Sin à + Cosg.Cos8. Cost 

jeśli wyrazimy czas w sekundach, a kąt godzinny w radianach, 
będzie: 


t = pt, gdzie u = ve , gdyż doba gwiazdowa 


zawiera 86164 sek. śr.; wstawiwszy (2) w wyrażenie (1), otrzyma- 
my po scałkowaniu: 


OC" 


(3) Q; = = © (Sine.Sinè.+, + Cosy.Cosô.Sint,) 


To Obliczymy z (2), podstawiając h = o, będzie mianowicie: 


(4) Cost, — — tg ẹ .tgò. 

Dla równika ę=0, y= Se zatem dzienna ilosé ciepta be- 
dzie w tym przypadku: 
(5) ES A Cos 8. 


Dla bieguna zaś, zważywszy, iż ę=—, To= ī: 


2Cu = 
(6) Qp = - Se Sin à. 


Dla natężenia promieniowania mamy wzór: 


E = = (Sin y. Sin à + Cos ẹ . Cos è. Cost) 


dla równika zatem: 


ER C Cos ò. Cos T 
dt |e=0 : 1 
a dla bieguna: 
do | C Sin à 
dt 9-7 LL 


wynika stąd ciekawy związek między natężeniem promieniowa- 
nia dla danej szerokości oe i promieniowaniem na równiku 
i biegunie: 

dQ | dQ 
dt BLIGE HA GC! A EWĘ e 


2 
T 


(7) 


Przez całkowanie dochodzimy stąd do związku między ilo- 
ściami ciepła: 


(8) Qo = Dan. Sin p + Q =- Cost 


Q—oznaczają odpowiednie ilości ciepła, otrzymane w tym 
samym przeciągu czasu. Oczywiście, iż taki związek nie będzie 
zachodził pomiędzy ilościami dziennemi ciepła. 

Posługując się wzorem (3), można obliczyć ilości dzienne 
ciepła dla różnych szerokości geograficznych i dla różnych poło- 
żeń słońca; następująca tabliczka przedstawia rezultaty tych ra- 


chunków, przyjęto tam Ce = 1000. 


Dla półkuli południowej, t. j. dla p < o, liczby tabl. I są 
słuszne pod warunkiem jednakże, aby 6 brać ze znakiem od- 
wrotnym. 

Tabliczka I pozwala dla każdego położenia słońca zna- 
leié odpowiednią ilość ciepła; jeśli położenie słońca jest wy- 
znaczone nie przez zboczenie, lecz przez długość ©, to à można 
obliczyć ze wzoru: 


(9) Sin © = Sine . Sin8, 
gdzie e jest nachyleniem ekliptyki (e = 2392731", 0). 


Tab. I. Tablica ilości dziennych ciepła na granicy atmosfery. 
Table de valeurs diurnes de la chaleur à la limite de l’atmosphère. 


Szeroko$ci geogr. — latitudes géographiques 
30° | 409 | 500 | 600 | 700 | 800 | 900 


m | 100 | 209 


+ | 9174 10147 1086,7| 1132,2| 1151,8 1149,4 1137,2 1175,1 1231,5| 1250,5 
+22 | 927,2 1017,6 1082,0| 1119,1| 129,7 1117,3| 1092,2 1105,9 1159,0| 1176,9 
+-20° || 939,7| 1020,6 1074,5| 1100,5| 109,0 107,3, 1032,0 1009,7| 1058,2 1074,5 
+169 | 961,8 102,1 1056,3| 1064,4| 1036,1 986,1 916,2 844,3, 852,6, 865,9 
+12 | 978,1] 1020.7 1033,6 10108 971,2, 899,2 805,4, 700,3, 643,2 653,2 
8 | 990,3 1013,6|1006,5| 969,7| 904,4 821,9] 699,2) 569,1] 445,3) 437,2 
40 | 997,6 1001,5 975,2] 919,4) 836,0, 727,4) 597,3 450,5 295,0 219,1 
1000,0 984,8 939,7) 866,0| 766,0 642,8] 500,0 242,0 173,6 0 

—4 | 997,6) 963,4 900,2 809,8) 695,1) 559,5) 407,5) 344,5) 79,1 
—8 | 990,3 937,6, 857,0) 751,1) 623,3, 478,0) 320,6 158,6) 14,7 
-—120 | 978,1| 907,3) 810,2) 690,2! 551,3 398,9] 239,8 86,4) 0 

—16 | 961,3 872,7 760,1! 627,4] 479,5 322,7 166,3 30,6 
—20° | 989,7] 834,0) 707,0) 563,2) 408,3 250,2 1015 0 

—22 | 927,2) 813,2 679,4) 530,7| 378,2) 215,7 73,0 0 

an | 917,4) 797,5 659,0| 506,9! 348,0 1915) 543 0 

| 


Zooczenia sloñca — déclinaisons dti solei 
=; 


SE (E => Te 


Aby ilości dzienne otrzymać w kaloryach, trzeba liczby po- 
EC się tyczy stałej ŻĘ jest 

1000 * SA 

to ilość ciepła, jaką odbierze w dniu porównania wiosennego 
lub jesiennego cm.? pow. ziemi na równiku, wynika to z (5). 
Przyjąwszy stałą słoneczną równą 2 kal. gr., zredukowaną do 
średniej odległości od słońca, mamy 


(10) log i = 2,96105 


wyzsze mnozyc przez 


Gdy chodzi o większą ścisłość, można uwzględnić zmianę 
7, mnożąc wartość znalezioną wyżej przez stosunek 


(|= | +e. ar 


r 1 — e? 


gdzie e — mimośród równy 0,0168, „ — długość perihelium. 


LN: AŻ 


Powyższy stosunek obliczać też można z szeregu: 


WER DI = 1,00014 -+ 0,03660 . Cos £ + 0,00071 . Cos 26 4. 
-+ 0,00002 . Cos 3$ +... 


$ — jest to anomalia średnia; $ = n (t — ty). 
Aby obliczyć ilość ciepła, otrzymaną w dłuższym przeciągu 
czasu przez cm.* na biegunie, wyjdźmy z wyrażenia: 


GQ. EE Gy. AE 
(12) KE Sin à — 78 Sine . SinO; 
Zmienna £ daje się łatwo wyrugować przy pomocy znanej 
równości z kinematyki dła ruchu środkowego (prawo pól): 
(13) r? d z =C 

dt 
c — jest to t. zw. stała pól. 
Otrzymujemy zatem: 


dOS EE ST 
(14) ROM Sin: . Sin O 
stąd, drogą całkowania, dochodzimy do ilości ciepła, odpowia- 
dającej okresowi czasu (©,, ©): 
(15) > A Sin e (Cos ©» — Cos O) 
O ma 


Jeśli szukamy ilości dziennej, to wyrażenie powyższe z po- 
wodu drobnej wtedy różnicy ©— ©, redukuje się do (6). 

Kładąc w (15) ©,—0, t. j. licząc ilość ciepła od porówna- 
nia wiosennego, będziemy mieli: 


(16) Qo=2£ Sin: : Sin ` 


Stąd dla ilości rocznej ciepła, otrzymanej na biegunie (w cią- 
gu !/, roku), wypada wzór: 


` C 
(17) SEA Sin e = 0,79620... -< ; 


jeśli chcemy wyrazić Q, w kaloryach, wypadnie wziąć da © 


wartosé, dana przez: 


C 523453 


(18) log -> 


otrzymujemy zatem w kaloryach gr.: 
(17) Qr, = 136626 kal. gr. (cm2); 


zwracamy uwagę, iż jest to ilość, otrzymana na granicy at- 
mosfery. 

Przejdźmy teraz do przypadku 9 = 0, t. j. do równika. 

Zadanie jest tu trudniejsze z powodu zmiany dnia i nocy. 
Ilość dzienną ciepła obliczać można ze wzoru (5). 

Co się zaś tyczy ilości, otrzymanej w dłuższym przeciągu 
czasu, postaramy się uczynić ją funkcyą ciągłą długości słońca 
i prostszą do obliczania, aniżeli metoda wskazana w rozdziale 
poprzednim. Weźmy w tym celu pod uwagę, zamiast ilości ciepła, 
którą otrzyma cm.?, ilość ciepła, którą odbierze pas szerokości 
1 cm., rozciągnięty wzdłuż całego obwodu równika. Oznaczmy 
tę ilość ciepła przez W; będzie to ciepło, otrzymane przez oświe- 
tloną połowę równika. Oczywiście W będzie zależało tylko od 
położenia słońca i będzie funkcyą ciągłą ©, gdyż ruch dzienny 
ziemi nie zmienia kształtu części oświetlonej, będzie on tylko po- 
wodował jednostajny rozkład ilości ciepła W wzdłuż całego pa- 
sa równikowego, którą to ilość w razie nie istnienia ruchu wi- 
rowego otrzymałaby tylko połowa równika. Wynika stąd, że je- 
śli W podzielimy przez obwód równika, to będziemy mieli śred- 
nią ilość ciepła, jaką otrzyma każdy cm? powierzchni ziemi na 
równiku w danym przeciągu czasu. Zrozumiałem jest, iż owa 
ilość średnia nie będzie ściśle odpowiadała rzeczywistości, za- 
wsze będą pewne drobne różnice dla dwóch miejsc na równiku; 
różnice te będą tem mniej znaczące, im większą ilość obrotów 
wykonała ziemia w danym przeciągu czasu. Teoretycznie rzecz 
biorąc, gdyby ziemia obracała się nieskończenie prędko, każdy 
cm.* otrzymałby jednakową ilość ciepła i równą ściśle ilości 
średniej, określonej powyżej. Następująca interpretacya wyjaśni 
najlepiej pojęcie średniej ilości ciepła. 


Przypuśćmy mianowicie, iż na powyższy pas równikowy pa- 
da pewna ciecz w takiej ilości, iż na każdy cm? pada tyle cm.*, 
ile kaloryi gr. odbiera dana jednostka powierzchni od słońca; 
przypuszczamy następnie, iż owa ciecz rozlewa się tylko wzdłuż 
pasa równikowego i dzieje się to nadzwyczaj szybko, tak że po- 
wierzchnia tej cieczy wszędzie będzie na jednakowej wysokości 
nad ziemią. łatwo zauważyć, iż owa wysokość nad ziemią bę- 
dzie właśnie równa średniej ilości ciepła, jaką otrzyma każdy 
cm? (jeśli jeden cm. odpowiada 1 kal. gr.). 


Aby obliczyć W, trzeba najpierw znać ilość ciepła dW, 
którą odbierze równik w czasie dł; powierzchnia ds odbiera 
ilość ciepła È Sinh. dt. de, należy więc wziąć sumę podobnych 


wyrazów i rozciągnąć ją na oświetloną część równika. Oznacz- 
my przez A łuk równika, zawarty między danym elementem do 
i miejscem na równiku, dla którego słońce przechodzi przez po- 
łudnik; widać, iż À będzie kątem godzinnym, mamy więc: 


f Sin h = cos à . cos À 


9) |  ds=Rdh 


R — promień ziemi w cm., A w radianach. 


Będzie zatem; 


a SH Cosh . dì; 
dt 

to: jest: ; 

(20) o = E cos 8; 


Jeśli przez dQ oznaczymy przyrost średniej ilości ciepła, 
czyli w naszej interpretacyi przyrost wysokości cieczy w cza- 
sie dź, to będzie: 


dQ 21) *44W 
(21) E Ed 


11 — 


dQ C 


zatem: Sady: cos à (22) 
Wyrugowawszy zmienną £, przy pomocy (13), mamy: 
dW C 
| ae = c 2R: Cosż 
(23) 
dQ "Lt 
| stars Cos à 


Zauważywszy, że cos 0 = V 1 — Sin’e. Sin? O ; 
otrzymamy średnią ilość ciepła w postaci całki eliptycznej nast.: 


IEC CS GT E 
(24) = SEN V 1— Sin’e. Sin © . do 
Oo 
Kładąc O = 0 i O= D mamy dla ilości rocznej wyrażenie: 
p SZARE: 
(25) Oz [*y1=Sinte. Sin ©. dO 


0 
Jest to znana całka eliptyczna zupełna 2-go rodzaju o mo- 
dule sine; oznaczają ją przez E (6). 
Wartości /og E (9), dla kąta modułu w odstępach Oft, po- 
daje prof. Witkowski w ,Tablicach Matematyczno - Fizycz- 
nych*. Mamy: 


Qr MA, = E (e) = > = E (23°2731,"0) 
po obliczeniu wypada: 
(26) Q: = 1,9182 . = 
w kaloryach zaś: 
(27) Q: = 329170 kal. gr. na cm? 


Jeśli teraz obliczymy Q;, sumując ilości dla poszczególnych 
dni (posługując się tabliczką w rozdz. I), to otrzymamy: 


Q; = 332640 kal. gr. 


RL E 


zaledwie o 1% różną od poprzedniej; różnicę tę można zresztą 
raczej przypisać niedokładnościom metody sumowań. 

Doszliśmy zatem do rezultatów, które dowodzą słuszności 
metody wyżej opisanej. Wartość (25) można przedstawić w po- 
staci szeregu, mamy mianowicie: 


RZA FREAK: KR: 
M 1—sin?e.sin O = 1 — Z sine. sin: © — — sinte. sin: O+... 


— À gäre, SON RE E - = 6. „m 3) inme, sinr© + 
Ponieważ zaś: AC 
(28) ER sint ©. dO = e 2 Se z: m — parzyste 
zatem: 
(29 Q = e [17 sine = Ę sint e — ZĘ: lte, rm 
ki Die mo E SE 


Postaramy sie teraz przedstawié wzôr (16) w postaci szere- 
gu Fouriera, któryby pozwalał dla każdego przeciągu czasu 
obliczać odpowiednią ilość ciepła; przyjmiemy ©, =o. Podob- 
nie, jak poprzednio, mamy rozwinięcie: 


GOA A PREZ A MARE 
dO ao! 5 Sin’ e. sin? O 


sin*e.sin* © —....) 


Gel — 


szereg ten przekształcimy, zmieniając potęgi wstaw przez do- 
stawy wielokrotności łuków; weźmy najpierw szereg ogólny: 


(380) S =a, + asin? x + asint x +... + am. sin™x +... 


mamy, jak wiadomo: 


Sin EPA 
Zä 


1 
Sin * x =Š — + cos 2x-+ + cos 4x 
8 Z 
à 5- 15 1 
H 6 = E Ee 
Sin *x 16 D cos 2x LA z COS 4X — zg COS 6x 


; 35 7 4 1 
Sin 3x= — À cos 2x + z cos 4x — cos 6x + — g COS 8x 
o 


16 16 
(31) 
63 49 15 157 
10y — "7 _ 
Sin!’ x = 256 256 cos 2x +67 cos 4x 515 COS OXT 
Peu cos 8x — 1 s 10x 
Gr ee EE 
231 99 55 
inl? y =. SS see, 
Sin!? x 1024 „zę COS 2x + gp COS 4x gjg COS 6x + 
> 8 - 10 12x; 
+ 1024 OR = 515 08 x+ ag” cos 12x ; 


Podstawiając te wyrażenia w (31), otrzymujemy szereg 
Fouriera: 
(32) S = A, FA; cos 2x + A, cos 4x + A; cos6x--.... 
współczynniki À, A, A... dane są przez szeregi następujące: 


A = atouts OR E Per ZEE 1 


8% 1167517128: 105679 "17" 
m! 
(2.4.6...m)*' 5 
Kia Lacka ESA Lo de RS 
1 RE Ae IE AT DORE 
33 1 owo masc EE 
( ) R. g 42 I 16 43 ‘ 39 4 l 54 *5 c.oo 
A SCH R 1 të 
z 82.70.16 br? 
A,= 1 a, ‘- 5 a. + 
128 Z EE 


Ee IAA" 


naturalnie, iż powyższe przekształcenie wymaga pewnych warun- 
ków co do zbieżności szeregu (31), którymi się zresztą tu 
nie będziemy zajmowali. Stosując powyższe wzory do rozwi- 


nięcia dQ , mamy: 


EN It (A, + À, cos 2O + A, cos 40 +....) 
c 


dO x 
A= ee 
z 4 64 166 cw 
1E 
F4 A, = 4 Sin? e + jg Sint «++ zy Eine. 
à l Sinte 8 one — 
Å = — gą Sin e 955 Sin Ne 


D 


Po podstawieniu danych liczbowych, wynika: 
(85) 45 = £ (0,30530 + 0,01818 Cos 20 — 0,00014 . Cos4© + ...\ 


Całkując, otrzymujemy szukaną ilość ciepła w postaci 
szeregu: 


(36) Q= 5 (0,30530 . © + 0,00656 . Sin 20 — 0,00003 . Sin 4O +...) 


Kładąc © = Oe, mamy ilość roczną zgodną z (26): 


C 
Q: = 1,91824. z 


Jeśli znamy przeciąg czasu, odpowiadający okresowi (0, ©), 
to długość ©, potrzebną we wzorze (36), otrzymujemy bądź z ta- 
blic, podających dla każdego dnia roku długość słońca i zbo- 
czenie, bądź też rozwiązujemy tak zwane zadanie Keplera. Szu- 
kamy przedewszystkiem anomalii średniej E: 


é=n(t—t); 


to —czas przejscia przez perigeum. 


= 


Se 1 B= 
Następnie obliczamy anomalię mimośrodową E ze wzoru: 
(37) E — e. sin E=$ 


stąd anomalia prawdziwa y otrzymuje się na mocy wzoru: 


(38) ez=y] Ee ez 


Jeśli w jest długością perigeum, to © = w+ Y. 
Można zresztą odrazu obliczać anomalię prawdziwą przy 
pomocy szeregu: 


(39) v=t+ pe $) Sin 1 Ze Sin2£+ äi Sin 3£ 4... 


podstawiając e = 0, 0168, mamy stąd: 


(40) ©=wv+t+-0,033599. Sin ¢-4 0,000853. Sin 2 £- 
4- 0,00005. Sin 3 £... 


Posługując się tym wzorem, możnaby też wyrazić Q w posta- 
ci szeregu: 


Q=A, +A6+A,Sin$4-B,. Cosé+A.. Sin26 +. 


pozwalającego bezpośrednio obliczać Q, znając odpowiedni prze- 
ciąg czasu !). 


Ill. Obliczenie ilości ciepła dla jakiejkolwiek szerokości geogra- 
ficznej 4. 


Obliczenie ilości dziennej ciepła podane było w rozdziale 
I-ym. Przystąpimy więc odrazu do obliczenia ilości średniej, 
postępując w ten sam sposób, jak poprzednio. Weźmy więc pod 
uwagę pas szerokości 1 cm., rozłożony wzdłuż całego obwodu 
danego równoleżnika. Jeśli, podobnie jak poprzednio, oznaczy- 
my przez X łuk równoleżnika, zawarty między punktem rozwa- 
żanym i punktem, dla którego słońce przechodzi przez południk, 
to będzie: 


1) Angot. Distribution de la chaleur sur la surface du globe. Annales 
du Bureau Central Météorol. Paris, Année 1883. 


LCE 


o f Sin À =Sin4. Sin à -|- cos #. cos à. cos À. 
| do—R. cosy. dA. 


Biorąc teraz całkę wzdłuż oświetlonej części równoleżnika, mamy: 


(2) He EE cos à. cos À) R. cosy. dh. 
0 
W — ma to samo znaczenie co poprzednio, X, jest połową oświe- 
tlonego łuku równoleżnika, daną przez wzór: 
(3) cos À, = — tg ọ. tgô; 
4, jest oczywiście równe połowie łuku dziennego słońca. Z (2) 
otrzymujemy : 


(4) joe RES „ Sing. Sin à. À, + cos ẹ . cos à. Sin à). 
Zadanie, jak widzimy, jest bardziej złożone od poprzednic- 
go, a to z powodu zmienności An, które dla równika było równe 


stale >. Zauwazymy, iż natężenie średnie promieniowania bę- 
dzie: 
GOL NZ AW 
dt ~ 2mR.cosy dt 
otrzymujemy: 
= (Sin ę. Sin 8.4, + cos%.cosó.SinA „) 
lub 


r Bi 
OR RĘ — 7 (Sinę . Sin 8.4,+- Cos ẹ . Cos ô . Sin à ] 
podstawiwszy E CE wypadnie po przekształceniu: 


(6) LS = EŃ = - (tg . Sin e. Sin O -^ + 


ść r 7 


Zadanie sprowadza się teraz do rozwinięcia wyrażenia po- 
przedniego. Drugi wyraz dwumianu w nawiasie rozwija się tak, 
jak w rozdziale poprzednim; aby rozwinąć pierwszy wyraz, trze- 
ba będzie znaleźć rozwinięcie A, Oznaczmy: 


S= Tgę. Sin e . Sin O.Ą,, ponieważ Cos À, = — tg ẹ . tg ò, mamy więc: 


= + arc Sin (tgę.tg 0). 


Stosując zaś rozwinięcie: 


H 1 3 > 5 5 NET] 35 3 
(7) arc Sin X =X + e * rio *t qm * +ga * + 
63 a 
| 2816" 
otrzymamy: 
z 1 3 
ào = 5 + (tg9.tg8) -+ (tee. tg) + qq (tg 9. tg 8) * +-.... 


założywszy dla skrócenia a =sine.sinQ), mamy też: 
Sine.SinQ 1 48) 5 


(00: m1 EE 2 g3 2 2 gó 5 11 
s /1— Sin? e. Sin © PON AŻ Ce 
WE SE 
ZE TIE SIN A 
Zatem: 


SE e PUR ON LE ee 
et eeh Ask Bebe Dee OBGĆ Ee 


21 89 1081 


| 1 3 3 3 5 gl 9 ja 11 
A 6 8 9. (a + o a + 4 A + 16 a Sé 128 a +.) + 


3 5 5 5 7 9 285 11 
+45 8 9 (8 ++ 97% + aa +...) ++ 


5 "ët 63 


Tep | 9 11, | 
- ro 8 F (a' -+ a+ at.) + 


„Prace“ Wydz. III. Ne 15. Badania teoretyczne ilości ciepła. 2 


a2 FES 


porządkując zaś według potęg a, będziemy mieli: 


(8) A = z+atgy+ać. E eer 18°? )+ 


3 l 2 
+o ZEE me) + 


Podstawiwszy to rozwinięcie w wyrażenie na S, mamy: 
(9) S=4Ąa.tgę=a, Sin © + äs Sin? O + az. Sint © + a,. Sinf © +... 


współczynniki a, 4, a;... są następujące: 
a, — 3 -tg9.Sine 


a= tg?o. Sin?e 


m 
w 
I 


(z wets tg? e Sint e 


EI 
a, = E tg? e+ n tot 9+ 70 w tg” G Lë: 
(10) \ 


a= [> Lee Lt +4 iere EIS 


35 
A= | iog t8 9+ A 9-318" 9 +3 z tę” pa 


ps Te tg" +) SE? 


z (9) wynika rozwinięcie Fouriera: 
(11) S=A,+A,.SinQ+-A,. Cos 2 O+ A,. Cos 4 O+-... 


A, =a, zaś współczynniki A, A, A, i t. d. otrzymamy, po- 
sługując się wzorami (33) rozdz. Il. Mamy też następujące roz- 
winięcie: 


Ee POJO ce 


3 Sie a. 
(12) y j (= Cosżę' Sin? O = B, + B, . Cos 2 © + B,. Cos 4 © +- 


+ B; . COS 6 © +... 
1Sinże 3 Sinfe 15 Sint: 
4 cos? 64costp 768 cosy `“ 
1 Sinż: 1 Sinte 15 Sinfe. 


Bt 


Be cosy | 16costę " 512 eoste | ™ 
(13) BSE Í Sint e 3 Sin$e 
4. 64costę 256 costę 
in 
SC 1 Sinfe 


= 512 costę | "" 


Wypada stąd rozwinięcie szukane: 


dQ cose C 
T 


EE 


(A, + Bo) + A; . Sin © + (A; + Bal Cos2© + 


+ (À, + B,) Cos 4 © + (A, + Be) . Cos6 O +- |; 


t. j. rozwinięcie kształtu : 


ANT 

(15) -==—(M,+M,.Sin O+ M, Cos 2 © + M,. Cos 4 © ++...) 
uw u 

godnem jest uwagi, iż dla szerokości —g, otrzymamy te same 

współczynniki M, M, M,..., za wyjątkiem M,, które zmieni znak 

na odwrotny, będzie mianowicie: 


(16) be EE EK 


doj 

Rozwinięcie to będzie słuszne też dla szerokości ę podczas 
miesięcy zimowych pod warunkiem, że © oznaczać będzie wtedy 
kąt liczony od punktu jesiennego «=. Jest to zrozumiałem 
też z tego powodu, że położenie słońca względem równoleżnika 
ẹ przy długości ©, jest takie same, jak dla szerokości —g, gdy 
długość wynosi t+ ©. 


VERT "WEN 


Jeśli liczyć będziemy ilość ciepła od punktu wiosennego, 
to z (15) wypada wzór dla Q: : 


C l 
(17) = < (M, +M, © — M, Cos© + 3. $in2 + 


L Sin 4 O+...) 

Oznaczmy przez Q+e ilość ciepła odebraną w okresie, gdy 
słońce jest w półkuli północnej, a przez Q—: gdy w południo- 
wej, będzie zatem: 

Cc, 
| Qt = (rM+2M) 
(18) 
C 
| Q—— 4 (1M,+2M,) 
Wypadnie też zauważyć, iż ilość ciepła, otrzymana w ćwiart- 


ce w=0,0= 2) równa jest ilości ciepła, odpowiadającej ćwiart- 


ce ës GES EX T ) to samo stosuje się i do pozostałych éwiar- 


tek. Dla rocznej ilości ciepła będziemy mieli wzór: 
dn M >> Cos p (A, 1 Bo) 


czyli: 


E 1 3 15 
(19) Q 2 Cose | CEE za rąk) 


(i 1 Sinře 3 Sinte 
TU 4 Cos? e 64Cost g 


Z (18) wynika również ciekawy związek: 


G, 
Q+:7 Q- = 74M, 


(EES 
C 
(20) Q++=Q—=2 —. Sing. Sine. 
Do wyników naszych zastosujemy dane liczbowe dla War- 


szawy; będzie zatem = 52,22. Wypadają następujące wartości 
współczynników: 


log a, =1,40122 ; a, = 0,2519 


| log a, = 1,90673 ; a, = 0,8067 
(21) | log a; = 2,51182 ; a, — 0,0325 


log a4 = 2,26002 : a, = 0,0182 
log a; = 3,62880 ; a, = 0,0043 


Na zasadzie tych wartości, otrzymujemy współczynniki: 


pp | Ao = 0,1449; A, = 0,8067; A, = — 0,1526 
(22) | À, — 0,0062; A, — — 0,0003 i t. d. 


Ze wzorów (13) otrzymujemy jeszcze: 


A B, — 0,8844; B, — 0,1187; 
Dis) += 0,0040; B, — 0,0001. 


Podsławiając we wzór (14) wartości powyższe, mamy roz- 
winięcie: 


(23) e G (0,2007 +- 0,1573 . Sin © — 0,0066 . Cos 2 © + 


+ 0,0004. Cos 4 O +...) 


Całkując, otrzymujemy szereg, przedstawiający średnią ilość 
ciepła dla Warszawy: 


(24) Q= Ze (0,1573 -+ 0,2007 . © — 0,1573 . Cos © … 


— 0.033 . Sin 2 © + 0,0001. Sin 4 O+...) 


wynika stąd: 
Q+-: = 0,9451 > Q—: = 0,3159. ©; 


zaś ilość roczna: 
Qr = 1,2610. =; 


W kaloryach kgr. będziemy mieli: 


25 Q+:= 162,20 kal. kg.; Q—e=54,21 kal. kg. 
| Qr =216,41 kal. kg. 

Zobaczymy teraz, iż szereg (23) wiąże się w prosty sposób z ilo- 
ścią dzienną ciepła, którą oznaczyliśmy przez Q;. 

Przy obliczaniu Qj za pomocą wzoru (3) rozdz. 2-go, przyj- 
mowaliśmy zboczenie 8 lub długość © słońca za stałe w ciągu 
dnia i równe wartości współrzędnych słońca w południe każde- 
go dnia. Q; będzie więc funkcyą ©; postaramy się rozwinąć 
ją na szereg. Tego rodzaju zagadnieniem zajmował sie Angot!); 
metoda przez niego używana, w celu odnalezienia rozwinięcia 
dla Q;, ma jednak charakter empiryczny. 

Zauważmy teraz, iż przy obliczaniu Qj mamy całkę iden- 
tyczną z całką, jaka występuje w wyrażeniu ilości średniej cie- 
pła. Jest mianowicie: 


BE ck (Sin  . Sin 8 + Cos ẹ . Cos à. Cos à) dh; 


Cosi =—tgy.tg0. 
(26) 


Q = ZP (Sin . Sinè-+Cosyw.Cosô.Cosr)dr; 
Cosr,=—tgy.tg0. 
wynika natychmiast stąd związek: 


(27) ` Uz Zen E R 


1) Angot (l. c). 


GE EE 


Rozwinięcie dla Qj otrzymać można zatem z rozwinięcia 
(23), mnożąc każdy -wyraz szeregu przez SECH Dzieki zatem 
tej okoliczności zadanie o rozwinięciu ilości średniej jest jedno- 
cześnie zadaniem o rozwinięciu ilości Q; traktowanem przez 
Angot'a. Należy jednak pamiętać, iż w szeregu, wyrażającym 
dQ 
dt 
znowuż w rozwinięciu ilości dziennej. Biorąc za r średnią od- 


ległość ziemi od słońca, mamy log ZE 3,72652; posługu- 


, długość © zmienia się w sposób ciągły, co nie zachodzi 


jąc się tą wartością, otrzymamy z (23) rozwinięcie ilości dzien- 
nej dla Warszawy: 


(28)  Q;=0,001 > (3,443 4- 2,698 . Sin O — 0,113. Cos 2 © + 
+ 0,007 . Cos 4 O --...). 


Możnaby też, na zasadzie rozwinięć (40) rozdz. 2-go, przed- 
stawić Q i Qj przez szeregi trygonometryczne anomalii średniej £. 

Obliczenia powyższe stosować można dla wszystkich równo- 
leżników, leżących poniżej koła biegunowego, t. j. dla szeroko- 
ści < 90°— = 66°33’. Dla miejscowości pasa biegunowego, 
rozumowanie poprzednie jest dopuszczalne . wtedy, gdy słońce 
oświetla tylko część równoleżnika, t. j. istnieje zmiana dnia i no- 
cy; dla okresu czasu zaś, gdy słońce nie zachodzi wcale, trzeba 
będzie rachunki przeprowadzić inaczej, niż poprzednio; będą one 
zresztą prostsze z powodu niezmienności À, = r. 

Rozpoczniemy od rozwinięcia ilości dziennej: 


(29)  Q= Ze, Sing. ginë 


r? 
godńem jest uwagi, iż Q;= Qp. Sins, Qp ilość dzienna dla bie- 
guna. Rozwinięcie Q; redukuje się w danym przypadku do je- 


dnego wyrazu, mamy mianowicie : 


Anm 
ŁNU 


7 Sing. Sin e. Sin © 


(30) El Ge 


OWA 2 


wzór ten stosować należy od chwili, gdy à=90°—% do chwili, 
gdy à, zmniejszając się, znów osiągnie tę same wartość, t. j. 
w granicach: 


DE Z OZGA RAC Ze 


A AT Sin e Sin e 


Według (27) będziemy mieli: 


10 sk Sin ọ .Sine. Sin O 


Jeśli oznaczymy przez Q ilość ciepła liczoną od chwili, gdy 


© 0»=arcSin Cos'ę 


Sin, | to będzie: 


(32) Que Sin 4. Sin e. (Cos O, — Cos ©); 


wypada stąd ilość ciepła, jaką otrzyma dany cm? w okresie cza- 
su, gdy słońce jest nad horyzontem: 


(33) OU — 2 Sin e. Sin e. Cosi, 


Oznaczmy teraz przez F (©) funkcyę, przedstawiającą śre- 
dnią ilość ciepła dla miejscowości pasa biegunowego, liczoną od 
©=0. Funkcya ta będzie się różnie wyrażała zależnie od po- 
łożenia słońca. Gdy © zmienia się od 0 do ©, istnieje zmiana 
dnia i nocy, jesteśmy zatem w warunkach zadania poprzedzają- 
cego; gdy © zmienia się od ©, do r —©),, funkcya schodzi się 
z funkcyą Q'+ q, gdzie q jest ciepłem otrzymanem w przedziale 
(0, Q,): dla © zawartego w granicach od t— ©, do t+ ©,» 
F (©) wyrazi się przez funkcyę, którą otrzymamy, dodając do 
funkcyi dla okresu (0, ©.) stałą Q's +q; wreszcie dla 


r-O<O<2r —O,, 


F (Q) jest stałą, bo słońce jest wtedy poniżej horyzontu i dla 
przedziału (2r—©,, 0) podobnie jak dla (0, CAL F(Q) przed- 
stawia się zatem, w 3-ch swoich fazach zmienności, przez różne 
szeregi; możliwe jest jednak znalezienie szeregu, któryby odpo- 
wiadał jej zmienności dla całego roku. 


dE 


IV. Straty promieniowania w atmosferze. 


Powyższe rachunki przeprowadzone były w tem założeniu, 
że promieniowanie słońca nie ulega zmniejszeniu przy przejściu 
przez warstwę atmosfery ziemskiej. Porównanie rezultatów otrzy- 
manych z doświadczalnymi wskazuje. na znaczny wpływ atmo- 
sfery na natężenie promieniowania i to tym większy, im słońce 
jest bliżej horyzotu, co się tłumaczy większą grubością warstwy 
absorbującej energię promieniowania. Rachunki, które mamy 
obecnie na celu, będą niejako drugiem przybliżeniem rachunków 
insolacyjnych; będziemy się tu starali uwzględnić tylko najwa- 
żniejsze czynności absorbujące, i ująć ich wpływ w możliwie 
prostą formę matematyczną. Dokładne uwzględnienie absorbcyi 
jest niezmiernie złożone i może być jedynie empirycznie trakto- 
wane. W powietrzu suchem głównym czynnikiem, który zmniej- 
sza natężenie promieniowania jest dyfuzya, czyli rozpraszanie się 
światła; obecność pary wodnej w powietrzu również znacznie 
osłabia promieniowanie. 

Oznaczmy przez Oe pierwotne natężenie promieniowania, 
przez q natężenie po przejściu drogi o w ośrodku jednorodnym; 
według prawa Bouguer'a natężenie będzie malało z postępem 
geometrycznym, t. j.: 


—ko 
(1) 47 9:€ 


k> 1 jest to pewna stała, zależna od ośrodka. Jeśli wzór (1) 
zastosujemy do natężenia promieniowania słońca, przypuszczając 
jednorodność powietrza i nie uwzględniając kulistości ziemi, bę- 
dziemy mieli: 

séi 

inh 

(2) FU 
h— wysokość słońca, p— współczynnik transmisyi, wyrażający, 
jaka część energii zostaje absorbowaną przy pionowem padaniu 
promieni. Wzór (2) wyraża tak zwane prawo Lamberta. Po- 
nieważ będziemy się nim posługiwali w następnych rachunkach, jako 
jednem z prostszych wyrażeń dla strat promieniowania, trzeba 
uzasadnić jego stosowalność do ośrodka niejednorodnego, jakim 
jest powietrze. Współczynnik k nie będzie wielkością stałą, lecz 


Zopa 


będzie się zmniejszał wraz z wzrostem wysokości H nad pozio- 
mem, a to z powodu malenia gęstości. Weźmy więc pod uwa- 
ge warstwę powietrza nieskończenie cienką dH, wzór (1) zasto- 
sowany do tej warstwy da nam: 


(3) dq =— K (H). q . ds 


lecz, jeśli przypuścimy, iż droga promieni jest prostolinjową, mo- 
żna przyjąć 


ó = ts do = RDA 
Sin h ? Sinh 
zatem: 
—_ KH) 
dq Sin h qan 
wynika stąd: 
H 


(MË rj a 
Log. (z) een Ten. 


U 


całkowanie należy rozciągnąć do granic atmosfery; q, będzie 
promieniowaniem na granicy atmosfery, równe, jak wiadomo, 
około 2 kal. gr., cm?, min. 

Jeśli założymy: 


H 
f xw „dH 
E 0 
(4) p=e 


wynika wzór Lambert'a: 


1 
Sin h 


(5) q=9.P 


zatem wzór (2) słuszny jest i dla ośrodka niejednorodnego, pud 
warunkiem, aby p brać ze wzoru (4). 

Wzór (5) był wyprowadzony .w przypuszczeniu, iż droga 
promieni słonecznych w atmosferze jest prostolinjową, nie uwzglę- 
dnialiśmy następnie kulistości ziemi; tego rodzaju przypuszcze- 
nie było dopuszczalne dla wysokości słońca blizkich 90°, dla 
wysokości jednak niewielkich pożądanem jest wyprowadzenie 
wzoru ściślejszego, w którymby powyższe dwie okoliczności były 


NEE 


uwzględnione. Weźmy więc pod uwagę element ds krzywej re- 
frakcyi i przyjmijmy, iż K=ap, v— gęstość powietrza, będzie 
zatem: 

H 


TA ` 
Lo EN Y ser ds 
P (z) k P 


0 
jeśli teraz złożymy: 


H H 
af ed fods 


(6) Kc NOZE 


foan 
o 


to będzie: 


(7) q=9.P ; (PS1; m=1) 


m będzie pewną funkcyą wysokości h, która dla h = 90° staje sie 

1; m jest długością drogi atmosferycznej. Wzór (7) 

jest wzorem ogólnym dla straty promieniowania, wzór Lam- 

bert'a jest tylko jego szczególnym przypadkiem, jeśli będziemy 
Se 


przyjmowac m = Gin 


Całe zadanie sprowadza się teraz do obliczenia funkcyi m. 
Oznaczmy przez i, kąt między styczną do krzywej refrakcyi i do 
danej warstwy powietrza. Niech 2 oznacza pozorną odległość 
zenitalną słońca, w średni współczynnik załamania dla danej war- 
stwy, R— promień ziemi; jak wiadomo: 

| (R+ H). Sin i = Const. = Ry,- Sin z 
(18) ie dH 
Cosi 


Wynika na zasadzie tego następujące wyrażenie dla drogi atmo- 
sferycznej: 


p dH ; 
(9 e |< — | Sin 27 
(3) d fa] Vi (R + ble "Sin? z 


RTE 


Z powyższego otrzymać można m w założeniu iż w—uw, czyli 
przyjmując drogę za prostolinjową. 

Posługując się rozwinięciami na szeregi, otrzymujemy. wtedy 
wyrażenie następujące: 


EX 
„/pHdH > nm 
S sec? Z. tg? z? 


R „/pdH R2 „fed 


podane przez Lamberta. Całki w wyrażeniu powyższem mo- 
żna obliczyć po przyjęciu pewnych hypotez co do zmiany gęsto- 
ści p powietrza z wysokością H. 

Przyjmijmy np., iż gęstość p maleje według postępu geo- 


(10) m =secz =śsecz. tg?z 


metrycznego t. j. p= Po - € Së , którą to hypotezę stosować mo- 
żna ściśle w razie równowagi izotermicznej; ponieważ 
H Be E 
8. fe dH = ciężarowi słupa pow. = p,, zatem: a= p, — L Lo 
0 


będzie wysokością jednorodnego słupa powietrza o gęstości py, 
który wywiera to samo ciśnienie p,. Wstawiwszy wyrażenie na p 
we wzór (10) i wykonawszy całkowania do granic atmosfery, bę- 
dziemy mieli: 

= 2; 1 Lo 2 La : 2 Lo tg” z 
(11) m=secz 9 R CZ: tg?z + (L, z RCos 2) 3R? Cossz 


Nie uwzględniając kulistości ziemi, otrzymaliśmy poprze- 


dnio m= SR =secz; widzimy z (11), iż błąd wten sposób 
1 
popełniony będzie mniejszy od = 1500 SCZ tg? z, gdyż L, = 10 km., 


np. dla h= 10°, biorąc m= -œ Sin $ Ti = ),/0, popełniamy błąd mniej- 


szy od — dla h= 45° błąd jest mniejszy od 205 , zatem nie- 
wielki. SE też w następnych rachunkach będziemy zawsze brali 


m = dla jakiejkolwiek wysokości, bo chociaż dla h< 10° 


sh nh 


błąd już nie jest zbyt mały, to jednak natężenie promieniowania 
jest wtedy na tyle słabe, iż przy obliczaniu całkowitej ilości 


KE TE 


dziennej ciepła błąd ten nie ma żadnego znaczenia. Dotąd, jakkolwiek 
braliśmy pod uwagę kulistość ziemi i niejednorodność powietrza, 
to jednak nie uwzględnialiśmy skutku tej niejednorodności t. j. 
skrzywienia drogi promieni słonecznych. Jednym z ciekawszych 
wyrażeń dla m, uwzględniających refrakcyę, jest wzór Laplace'a, 
którego wyprowadzeniem zajmiemy się obecnie. Weźmy pod 
. uwagę łuk MM drogi krzywolinjowej promienia; z rysunku 
widać iż: 
M 
Refr. = | du, gdzie u oznacza kąt, 


Mo 


który tworzy styczna z prostą prostopadłą w punkcie M,. Na 
mocy związków: 
s do AH e G 
| tgi=r-—; ruSini=R.p, Sinz; du = di + dt 
(12) ui 
| r=R+H P 


} : i w DBA, OC ER 
otrzymujemy : ei + w Cotgi.di= 0 [R RAR 24 


par rr, g 
Stad: SE _ Cotgi. du; i Fi 
Rn, R ' d 
Tu d D 4 
du=— R v„|2 x m $ 
y` — K le] „Sinżz | 7 
mm Kaes 
Wynika stad podstawowy CN 
wzór refrakcyi: z 
0 * 
f R Ho Fig. 1. 
j J e dy. 


(13) Refr. = = är 
J y 1- (sv | Sintz t 


Granicą RES jest 1, By dla SA =l 


Dr EE 


Przyjmijmy teraz według Laplace'a: 
m =lF2cp 
H 
c— pewna stała, mamy pdp.=cdp, pdH=—L,dp; /pdH=p,L,; 
o 


podstawiwszy te wartości we wzór (9) otrzymujemy: 


WA 
w du 


1 
Ke | y -- (Fe) ae 
si ry 


z porównania (14) z (13) wynika: 


_ 1 d(Refr.) , ry 
(15) dm cp, Sinz BR 


Przyjmując teraz u? Ry; = 1 (w rzeczywistości bardzo bliz- 
o 


kie jedności), otrzymujemy wzór Laplace'a 1). 

E ME, 
(16) MK. Cosh’ 
; macl 
Dla statej K mamy wartość ee 
jeśli Refr. wyrażamy w sekundach, wypadnie przyjąć: 

K = 58”, 36. 
Np. dla h= 10°, Refr.=320” mamy m = 5,63, tego samego 

rzędu co wartość przybliżona m = zz =5,76. 


Wszystkie poprzednie wyrażenia wyprowadzone są w przy- 
puszczeniu niezmienności temperatury, zresztą wpływ tej osta- 
tniej, jak również i wpływ refrakcyi jest w naszych rachunkach 
nieznaczący. Bemporad?) ułożył tablice wartości m dla 
każdej wysokości, biorąc pod uwagę zmiany temperatury i do- 
kładniej uwzględniając refrakcyę promieni. 


1) P. S. Laplace: Mécanique Celeste IV. Chap. 3. Paris 1805. 
*) Bemporad. Zur Theorie der Extinktion des Lichtes in der Erd- 
atmosphäre. (Mitteilungen der Grossherz. Sternw. zu Heidelberg). 


eet Bf 
Co się tyczy prawa Bouguer’a, wyrażonego we wzorze: 


q=q,.p lub logq=a — bm 


nie jest ono ścisłe w zastosowaniu obecnem. Bemporad!) 
proponuje przyjąć: 


a 
log q= a—bm 


gdzie współczynniki wyznaczyć należy doświadczalnie. 

Bartoli podaje też następujący wzór empiryczny: 

—a 
4 q. M 

Na zakończenie naszych rozważań ogólnych wypadnie nad- 
mienić, iż miały one charakter mało precyzyjny, a to choćby dla- 
tego, iż samo pochłanianie braliśmy zbyt ogólnie, nie wdając się 
w analizę oddzielnych części widma, dla których przepuszczalność 
powietrza jest niejednakowa; przyjęcie więc stałego współczynnika 
transmisyi p nie odpowiada rzeczywistości; rachunki oparte na tem 
przypuszczeniu będą więc tylko przybliżeniem rachunków inso- 
lacyjnych, których większą lub mniejszą dokładność będziemy 
mogli ocenić w ostatnim rozdziale. 


V. Obliczenie ilości ciepta dla równika z uwzględnieniem strat 
promieniowania. 


Opierając się na prawie Lambert'a, możemy napisać dla 
natężenia promieniowania następujące wyrażenie: 
1 
dQ; C A Sin h 
(1) ER 3 = ER 


Zauważywszy, iż dla równika sin h = Cosà.Cos-, otrzymujemy 
dla ilości dziennej wyrażenie: 
Fi ee WE 
f! Cos å. Cos T 
Cnu. 
(2) Q = Cas crus de 
] 
“o 


1) Bemporad. L’Assorbimento selettivo della radiazione solare nell’ 
atmosfera terrestre (Reale Accademia dei Lincei). Roma 1908. 


Ee E 


zaś dla natężenia promieniowania średniego: 


C D Cos à. Cos À 
Ss | Cos à. Cos A.p . dh 
] 
6 
A ma to samo znaczenie co poprzednio. 

Całki powyższe można obliczać jedynie metodami przybli- 
żonemi. W ten sposób postępował Angot (l. c.); otrzymawszy 
szereg wartości dla Q; dla różnych pór roku, szukał następnie 
rozwinięcia Fouriera metodą empiryczną, co zresztą dało re- 
zultaty zadawalające. Naszym celem będzie dojść do rozwinię- 


cia Qj lub ie przy pomocy metod analitycznych. Wprowadza- 
jąc do (3) zmienną ©, mamy: 

SecA 
V 1— Sin? s. Sin © 


d Gë ees Ae 
(4) A LE | V 1 — Sin? e Sin? O.p 
0 


„Cosa. da 


Cała trudność zadania polega teraz, jak widzimy, na rozwinięciu 
funkcyi 


Sec À 
V 1—Sin’e. Sin © 


P 


na szereg. Załóżmy: 

(5) x=Sine.SinO; «a = Seca. logo > 

Mamy zatem rozwinąć na szereg Taylora funkcyę: 
SES 

(6) E ; a>0. 


Dla x=0 funkcya ta, jak również i wszystkie jej pocho- 
dne, jest ciągłą; warunki konieczne rozwijalności na szereg Tay- 


lora są więc spełnione. Punkty x=+ 1 są dla niej punktami 
krytycznymi; wewnątrz koła o promieniu 1 niema żadnych pun- 
któw szczególnych, a ponieważ funkcya (6) jest parzystą, winien 
więc istnieć szereg kształtu: 


a 


HAE 


(CZY re = ag i aj x EAN TRAĆ Raz r; | ZJ 


zbieżny dla |x|<1 i któryby w kole danym przedstawił funkcyę 
(6). Aby dowieść iż szereg (7) będzie przedstawiał daną funkcye, 
należałoby najpierw wyszukać wyrażenia dla reszty szeregu 
Taylora, odpowiadającego funkcyi (6) i następnie dowieść, iż 
reszta ta dąży do zera. Przy sposobności wypadnie nadmienić, iż 
sam fakt zbieżności szeregu (7) nie dowodziłby jeszcze, iż gra- 
nicą tego szeregu będzie funkcya (6), wystarczy choćby wziąć 
1 


pod uwagę funkcyę f(x)=Sin x+ e * , której rozwinięcie, jak 
łatwo się przekonać, będzie identyczne z rozwinięciem sin x, bę- 
dzie zatem zbieżne, a nie będzie jednak jego granicą f(x) lecz 
sinx. Granicę dla reszty szeregu (7) znajdziemy przy pomocy 
teoryi funkcyj zmiennej zespolonej. Mamy więc funkcyę: 


rz: Vi — z’ 
(8) f(z)=e 


holomorficzna dla |z|<l. Całkując następnie obie strony toż- 
samosci: 


f(z t(z) f (z) f (z) WK (z) m f (z) 
2 X y + 2? E (Lee) X ZET, SEA zelt x) 


wzdłuż koła C o promieniu jedność (z— punkt bieżący na okrę- 
gu, x — punkt stały wewnątrz koła), mamy: 


F(x) = f (0) + zf (0) +-....-+ zj PO CAB ES 


zn+1 (z—x) 


„Prace“ Wydz. III. M 15. Badania teoretyczne ilości ciepłą. 
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Mamy zatem wyrażenie dla reszty: 


3 d xn+l / f(z).dz 
(3) Dn äs"? J za+1 (zx) 


Dowiedziemy teraz, iż Rn zmierza do zera, gdy m dąży do 
nieskończoności. Wiadomo, iż moduł całki funkcyi danej, wzię- 
tej wzdłuż krzywej C, jest mniejszy od iloczynu z długości 
luku C przez największą wartość modułu funkcyi danej na łuku 
C, będzie zatem: 


t | e (z).dz 


Sin? f(z) 
(10) | Ra | — "Oe | J BCEE 


= |x|n+1, Max 
kl zat! Cal 


Jeśli zauważymy, iż |z|=1; |z—x|>1— |x|, mamy: 
(11) RASE E Max |{(2) | 
Ge Che 
Należy obliczyć moduł funkcyi f(z), założywszy z=é+2ln,tin 
współrzędne punktów koła zbieżności, otrzymamy: 


a 


Papi — „760 H 


f(z) — e Vi = e 
stad: 
äu ER 
|i(z)|=e Wąs Ve a 
a i 1 
AP z ui ni U 
(12) |(z)|=e 
zatem, dla q — 1, mamy: 
z zin TT NTS 

(13) Maxje VIP? |=e ` 


otrzymujemy stąd następującą granicę dla reszty: 


x|n+1 «|/2+2V2 
Ra |< SC? 


(14) qe 


9585 Ee 


ponieważ |x|<|, oczywistem więc jest, iż dla wszystkich punk- 
tów wewnątrz koła zbieżności będzie: 


lim | Ra = 0 
Dowiedliśmy zatem, iż szereg (7) jest zbieżny i jego gra- 
a 
KS y 1—x 
Granica (14) dla reszty posiada dla nas ważne znaczenie 
praktyczne, pozwala bowiem ocenić wielkość błędu popełnionego. 
Przejdźmy teraz do obliczania współczynników a, szere- 
f@) (0) 


gu (7); możnaby posługiwać się wzorem A sposób ten 


nicą będzie e 


jednak nastręcza wiele trudności natury praktycznej w obliczaniu 
pochodnych; użyjemy tu innego prostszego sposobu. Dla uła- 
twienia będziemy szukali współczynników rozwinięcia: 

Ki a 


Vi x 


Vi—x e : H 3 
e , które będą te same co i dla e , mamy zatem: 


a 


paz iż à Ri E Eege 


dy CEVI RAS pa Vil ACC = 

stąd. Ska e Gun en zaju% 

=- la, Fax+ ax +. (rat DE ap A oj 
Daję Was X Š ; 128 


693 ,. 3008 „, , 45045 +; 765765 s, | 
"956" "1024" | 14336 * | 2298767 


Jeśli pomnożymy dwa ostatnie szeregi i przyrównamy współ- 
czynnik przy jednakowych potęgach zmiennej x ze współczyn- 
nikami rozwinięcia: 


ay = a, + 2a,x + 3a3x? + 4a,x* +-..... 


to bedziemy mieli: 


E 


—a 1 l S 
de ; 4 er a ao; a, = — paly aota] 
L (8 3 
a; — 5% a [ga atga ta) 
En ga DEE -a F- yata) 


EE, SE AWZ 
(ONES lt 198% ' 16% e para) 


pe (898249305 edo ND SA A 
% 137 {56 % 128% | 16% 78 % 19% d 
3-5. (8008, 7,7608 4 m8160% 0852016 
EE Lon" 256% KE Ee 


2 
"ais 


Jeśli teraz kolejno, poczynając od ay, wstawimy wartości 
współczynników w następne, to otrzymamy wyrażenia dla an 
w postaci wielomianów stopnia m—1 względem o, będzie mia- 
nowicie: 


ao = e, 
|| 
a= geer 
a=- (8—0) 
ae 2 
as =——48 (15 — 9a + a?) 
W a,=— Ser (105 — 87a. + 18a? — 08) 


kap W L JE 2—30a3 + at 
a; 3840 (945 — 975a +- 285a? — 30a? + at) 


(23562 — 37935a + 140400? — 211543 + 
+ 185at — 3a5) 


FEE EOM 


EE AE 
Stosujac otrzymane rozwiniecie do (4), mamy, jesli: 


x=sine.sin©; a= sech. loga = 


Hu 
zc dQ I fa 
C eeng el sin ©. d a, . COS À . dA + 
8 8 g 
sin? e. sin ©. | a cosà. dA + sine, sin* O. E a, . cos À. dA + 
s 0 


3 Ò 
Ep 
da RA mi da.) 
O 


rozwijając pierwiastek i mnożąc szeregi, otrzymamy ostatecznie 
szereg: 


(18) O =4, +, Sin O+ . sin* O-+%,. sin*O+..... 


gdzie założono: 


Koh 
y=- sin? y (b +1) 
ją = - (I, + 1, + I) 
(19) 
Eee gg Gl +121—155—1,)- 
y= 2 Zus, + 151, — 451, + 141, — 1) 


Wielkości I, l, L.... są dane przez całki następujące: 
ER 
= | "sch . cos „ dh; I, = loge 


n 
| ” psech ? da 
o o 


o | = 


(20) 


L = loge | “sec À. ps. dh...; kg "secm-1Ą, psech dh; 
o o 


EE 


Z (18), dochodzimy znanym sposobem do rozwiniecia Fourier’a: 
ANSE E 

(21) ay =— z (Ao + A, . Cos 20-+ A; . Cos 4O +....) 

współczynniki mają wartości następujące: 


EE 
At gh! ga | gt "ost 


A — se 
(22) "=" 9 Tu 9 Ta 39 Ta 16 14 
1 3 7 
= = D E — DW CH 
A 8 12 16 13 H 39 ‘4 .... 


otrzymane ze wzoru (33) rozdziału 2-go. 
Z szeregu (21) otrzymujemy natychmiast ilość dzienną cie- 
pła, podobnie bowiem jak w rozdziałach poprzednich mamy: 


=. r? d ) 


8 86164 
c — stała pól; bas > zatem: 


E 


(23) Q = A + A, . Cos 20-+A,. Cos4©+....) 
Całkując teraz szereg (21) dochodzimy do wyrażenia średniej 
ilości ciepła, jaką odbierze cm? powierzchni równika w okresie 
czasu (0, O): 


(24) == (Ao „O+ Sin20+ Žž, Sin4O+.. > 


należy pamiętać, iż wzór powyższy wyprowadzony jest w zało- 
żeniu, iż współczynnik p jest stały wzdłuż całego równika (w rze- 
czywistości waha się między 0,7 i 0,8). 

Z (24) 0 ilość roczna: 


3 5 |, 35 221 
(25) Q= n+- ht ghtjgit ue LT 


ee 


Aby sprawdzić rezultaty otrzymane, szukajmy na co się za- 
mieni wyrażenie (25) w przypadku p — 1, t. j. w przypadku roz- 
ważanym w rozdziale 2-im, gdy całkowita ilość ciepła dochodzi 
do danej jednostki powierzchni. 

Przedewszystkiem widzimy, iż dla p = 1, wyrażenia /m sta- 
ją się nieoznaczone: 


3 i” sec™—1) Dësch, dì] —=0. 
oli 


o p—1 


E 
(Im pr = (96, p 


gdyż | secm-14, dA = co. Zadanie sprowadza się więc do zna- 


o 
lezienia wartości funkcyi: 


A | 
secm—lx, zseA, dA 


Fu (= log, z fie sec™—1) . zsee à, dÀ ES =— 
A SE S (loge z)-" e 


O 
dla z—1. Według znanej zasady Æ(1) będzie równe stosun- 
kowi pochodnych, dla z — 1, licznika i mianownika, zatem: 


d f secn-th. zeeen „da G 1 | "semi. zsec EU 
dE dxo z. d zo 


d -m ; -m—1 l 
dz (loge z) l, | | SH log Z., 28 


dej 1 SA 
Fm (1) mm, Zi Fmi 0) 
Wynika stąd, iż F.=F,=F,=....=0, 
z wyjątkiem: F, =f? così. dì = 1; 


o 
zatem: L=1; L=k=1l, =.,=0; 


Będzie zatem według (19): 


Wel UE } sin? e; % = — -z sinte it. d. 


CAES 
stad w przypadku p = 1 otrzymujemy: 


_2C LE gr Eé x 
de RFO GAN STOSS AR 
szereg ten jest identyczny z szeregiem (29), otrzymanym w roz- 
dziale 2-im; potwierdza to słuszność szeregu (25). 

Przejdźmy teraz do danych liczbowych. Przyjmijmy np. 
p=0,8. Trzeba będzie najpierw obliczyć wartości catek I, I, l... 
danych przez (25); można tutaj użyć bądź metody graficznej, 
bądź też stosować rozwinięcie (7). Będziemy mieli np. dla I, 
szereg następujący (SInA= pl: 


1 1 
1 , 
h= J DH WE atiati a, ET ovo 
ag Ay As .... Są to współczynniki rozwinięcia (7), otrzymane ze 


wzorów (16) lub (17), jeżeli w nich podstawić p= 0,8. Be- 
dziemy mieli: ay, = p = 0,8; a, = — 0,0892; a = — 0,06196; 
a; = — 0,0485; a, = — 0,0403; a; = — 0,0345; a, = — 0,0203; 
it d. Stąd wypada: I, = 0,7270. 

Dla I, i I, wypada z rachunku graficznego: 


l, = 0,2204; I, = 0,0825; 


poprzestaniemy w naszych rachunkach na tych 3 wartościach, 
dokładność bowiem będzie wystarczająca. Wstawiając otrzyma- 
ne wartości we wzory (19), będziemy mieli: 


Yo = 0,7270; yı = — 0,0730; Ya = — 0,0031; 
Stąd na zasadzie (22): 
A, = 0,6894; A, = 0,0380; À, = — 0,0004; 


zatem, dla pochodnej średniej ilości ciepła, w przypuszczeniu 
p = 0,8, otrzymujemy następujące wyrażenie: 

AG 1 o 
(26) dO" z e (06894 + 0,0380. Cos 20 — 0,0004. Cos 40 +... 


RAT EE 


stąd, według (23), ilość dzienna ciepła będzie: 


90. 
(27) Q;=  ,"(0,6894 + 0,0380 . Cos 2 © —0,0004. Cos 4 O + ...). 


Jak już wspominaliśmy, Angot (1. c.) dochodzi do wyrażeń ilo- 
ści dziennych przez szeregi przy pomocy metod interpolacyjnych, 
otrzymuje on w ten sposób dla p=0,8, szereg tego samego 
kształtu: Q;= 689,6 + 39,1 . Cos 2 © — 0,4. Cosic, gdzie przyjęto 
2Ch 


T? 


= 1000; widzimy zatem, iż różni sie on niewiele od (27); 


drobne różnice współczynników należy przypisać niedokładno- 
ściom rachunków. 
Z (26) otrzymujemy wreszcie dla ilości średniej szereg: 


Q= 5 (0,21944 . © + 0,00604 . Sin 2 © — 0,00003 . Sin 4 © — ...) 


kładąc © = 2x, mamy ilość ciepła, otrzymaną na równiku w cią- 


gu roku: Q, = 1,3788 . £ = 236,610 kal. kg. cm?. Angot otrzy- 


muje metodą sumowań 236,1 kal. kg. 


6. Obliczenia dla bieguna z uwzględnieniem strat promieniowania. 


Wyrażenie dla ilości dziennej otrzymujemy . natychmiast 
z uwagi, iż na biegunie wysokość słońca w ciągu jednego dnia, 
trwającego 24 godziny, można przyjąć za stałą i równą h= ð, 
będzie zatem: 


1 


Si 
(1) Or es E Sinó.p Ce à — zboczenie słońca 
lub tez 

C . C 
un Q= HEI SineSinO.p D 


Ruch dzienny słońca niema żadnego wpływu na ustonecz- 
nienie na biegunie, tak iż nie trzeba wprowadzać pojęcia ilości 
: 2 


Mala? ` am 


średniej ciepła; ilość ciepła, odpowiadająca okresowi (O, O), 
będzie daną przez całkę: 
223 
o Q=% (Sine.SinO.p 
o 


Cosc e . Cosc © 
AS 


zaś ilość roczna: 


Q: de [?Sine.Sin©.p 


o 


Cosc e . Cosc ©) 
. „dO; 


gdyż biegun północny (resp. południowy) oświetlony jest w okre- 
sie (O, x) resp. (x, 2x). 

Przy pomocy przybliżonych kwadratur, można będzie z (2) 
obliczać Q dla każdego ©. Postaramy się jednakże, podobnie 
jak poprzednio, podać Q w postaci szeregu Fouriera; trzeba 
będzie rozwinąć: 

1 a 


Sin e Sin. © "VI geet e 
p =e V OS (©) Ee loge 


Według rozdziału poprzedniego będzie: 
Vi Coeie 
(4) e = à, + a, . Cos? © + a, Cost © + a; . Cost O +... 
e / sł R dO : 
Zanim wyprowadzimy rozwiniecie dla iS: zauważmy, iż 


z (4) otrzymać można zaraz rozwinięcie ilości rocznej, jest mia- 
nowicie: 
T 

2C LAVER 

CO. = = Sine f (a + a, Cos? © + a, Cost O +...) d (Cos O) = 
o 
1 
= sine f (aa + a, X? + 4 x* +...) dx; 


O 


Zatem 


6) Vo Be gi...) 


wyraża sie więc w prosty sposób przez współczynniki rozwinie- 
cia (4). 

Z (4) wynika, po zastąpieniu potęg przez wielokrotności 
łuku: 


A vi — Cos? © 
(6) e = Be + B; . Cos2 O + b, . Cos4O + 8;.Cos 6 O+... 


Bos Bu Bae otrzymamy ze wzorów rozdziału 2-go, zastępując x przez 


SX Mnożąc (6) przez © Sins, Sin © i słosując tożsamość: 


Sin. Cos m O= > Sin(m+1)0- + — DG 


otrzymujemy wyrażenie dla pochodnej: 


c 
(7) a Sine | (6 — EU TE Fa Sin 3 O + 


Ab Sin5©+ „| 
po scałkowaniu zaś: 
C i K Ba B e B ś 
Q= € Sin: | (6 — 3 5-57 .|-- (6 lee? 
„ko -fa Cos3 ©— bo ba Cos5O+.. |: 


t. j. kształtu: 


8)  Q=--Sin:(A,+A, Cos®+ A, . Cos 3O+ A, Cos 50+...) 


da Lee 


Współczynniki A, A, A... mają wartości następujące: 


Ri = us EA 5 a E 

1 EE 

Ee EE Ee E E 
(9) Aa” Bag 92 64 e q92 Bt ee 0 

n=- 1 a, — a ZI pa az" 

80% 64% 64 * 512075 
A SZR y — l REPAUA: = 

17 448 *8 256 *4 1024% 


Łatwo sprawdzić iż otrzymana stąd ilość roczna 
Q, = = SineA, 


będzie zgodną z (5). Zastosujmy teraz dane liczbowe; przyjmu- 
jemy, jak poprzednio, p=0,8; będzie więc 
1 


g = = 
Sine 


loge a= 0,5605. 

P 
Postugujac sie wzorami (16) rozdz. poprzedniego, otrzymamy 
kolejno: 
log (+a) = 1,75656; a, = 0,57090 
log (— a) = 1,20414; a, =— 0,16001 
log (— a) = 2,98937; a, = — 0,09758 
log (— a) = 2,83548; a; = — 0,06846 
log (—a,)= 2,71122; a, =— 0,05143 
log (—a;) = 2,60463; a, = — 0,04023 
log (—a,) 2,52023; a, =-— 0,03313 
log (— a) = 2,34209; a, = — 0,02198 
log (— a,)= 2,82689; a, = — 0,02122 
log (—a,) = 2,22667; a,=— 0,01685 


(10) 


GARE 


Współczynniki a, maleją bardzo wolno, tak iż w następ- 
nych rachunkach trzeba będzie użyć większą ich liczbę aniżeli 
w (10). Aby sprawdzić wartości (10), zastosujmy rozwinięcie dla 


a 


ei WE a 
e 


przy x= V3” bedzie zatem: 


1/1 xs d 
Ken EC 
$ TRE Se RR BEST WNE) EE EE 


5 V 6 
podstawiajac (10), otrzymujemy e = 0,503327 . 
Aby obliczyć granicę błędu, zauważmy iż składa się on z dwuch 
części: błędu, wynikającego z odrzucenia części dziesiętnych i na- 
stępnie błędu podanego przez granicę reszty 


|xjn+! -5V 242V 2 
(a? | 


RA < 


w danym przypadku n=22, gdyż wyraz 11-ty nie wpływa na 
szóstą cyfrę, zatem "RA < 0,000003, błąd zaś całkowity będzie 
< 0,000013; rzeczywiście, wartość dokładnie obliczona jest 


— a «/ 6 
e 2” =0,503325, 
różni się więc tylko 5-ą cyfrą od wartości powyższej. 
Mając wartości (10), otrzymamy ostatecznie rozwinięcie dla 
bieguna: 


dl) Q= = (0,18691 — 0,20117 . Cos © + 0,01137 . Cos3 O+ 


+ 0,00221 . Cos 5 © + 0,00039 . Cos 7 © + 0,00020 . Cos 9 © + 
-- 0,00007 . Cos 11 © + 0,00002 . Cos 13 O +... ). 


Gel 


Dla bieguna południowego rozwinięcie będzie to samo, 
pod warunkiem, iż za © należy brać © —#x. Podstawiając w (11) 
= z, otrzymujemy ilość roczną ciepła. 
C 


c 


(12) Q = —. 0,37382 = 64,150 . kal. kg. 1 cm?. 


Rozwinięcie ilości dziennej ciepła: 
ES 
P d do 


bedzie zatem: 


OC 


7 (0,632 . Sin © — 0,1072 .Sin3 © — 0,0347 . Sin5© + 


(13) Q= 


r 


— 0,0136 . Sin 7 © — 0,0057 . Sin 9 © —...). 


7. Obliczenia dla jakiejkolwiek szerokości 4 z uwzględnieniem 
strat promieniowania. 


W przypadku rozważanym, stosując wzór (1) rozdziału (5), 
mamy dla ilości dziennej wyrażenie; 


To 
(0) Q= A | (Sin q . Sin è+ Cos ę Cos ê Cos). 
o 


1 
Sin p. Sino-|- Cos ẹ . Cos è . Cost 
„dr 


P 
zaś dla pochodnej ilości średniej ciepła: 


4 


EC 
(2) zż=— © J Sing. Sin8-- Cosg, Cosa Cos). 
O 


1 
Sin o . Sint Cos g . Cos 6. Cosi 
p „dh 


Cos XA, = Cost 7—tgyę.tgó. 


TE 


Z (2) wynika ilość ciepła, odpowiadająca okresowi (O, ©): 


O b 
(e . m 

EGE << | | (Sing. Sinè+ Cosę. Cos à. Cos). 
ò 6 


1 
Sin e Sin 8-|- Cos ọ . Cos à. Cos À 
p d. dO. 


Podobnie jak poprzednio będziemy się starali rozwinąć (2) 
na szereg Fouriera. Postępowanie jednak w taki sposób, jak 
poprzednio, t. j. za pośrednictwem szeregu Taylora, chociaż 
teoretycznie możliwe, nastręcza jednak prawie że niepokonalne 
trudności natury praktycznej; trzeba będzie więc użyć innego 
sposobu. Angot, w podobnych przypadkach, używa metod in- 
terpolacyi peryodycznej, co daje zresztą rezultaty zadawalające; 
my postąpimy tu drogą odmienną, oprzemy się mianowicie na 
teoryi szeregów trygonometrycznych. Wiadomo, iż każda fun- 
kcya może być w danym przedziale (a, b) przedstawioną przez 
szereg trygonometryczny, pod warunkiem, aby w tym przedziale 
posiadała punkty nieciągłości pierwszego rodzaju; są to warunki 
wystarczające (lecz nie konieczne), znane pod nazwą warunków 
Dirichleta. Jeżeli zatem mamy funkcyę F(x), to w przedzia- 
le (a, b), może być ona przedstawioną przez szereg trygonome- 


tryczny łuku t= (2x—a— b) (przedziałowi a, b zmiennej x 
odpowiada przedział —x, x zmiennej t): 
4) Fœ@= p= ` + (a, Cost +b, Sin t)+ (a, Cos 2t+ 


+ b, Sin 2t) +... 


współczynniki a, a, Dän DA i t. d. dane są przez wzory: 
Kä CH 
Onea ee © (t).dtzam=— | Cos mt. ®. dt; 
j J 


ee 
Dm — = | Sin mt. gp (t). dt; 


2 MADA CH 


Nie wdając się w szczegóły, nadmienimy jeszcze, iż rozwinięcie 
(4), przedstawiałoby funkcyę F(x) na całym obszarze jej zmien- 
ności tylko wtedy, gdyby F(x) była funkcyą peryodyczną o pe- 
ryodzie równym (b— a). 

dQ 


Zastosujmy powyższe wzory do rozwinięcia funkcyi d 


dQ 


o peryodzie równym 2m. Funkcya AS posiada tę własność, iż 


[Oh 


przybiera wartości równe dla położeń słońca symetrycznych 
względem linji „rak—koziorożec*; wyraża się to jak następuje: 


e Pä. 


łatwo dowieść na zasadzie tego, iż współczynniki asm41 i Dom 
t.j. współczynniki wstaw ze wskaźnikami nieparzystymi i współ- 
czynniki wstaw ze wskaźnikami parzystymi będą równe zeru; 
oprócz tego, ze względu na symetryę, wyrażoną przez (6), cał- 


kowanie wystarczy rozciągnąć na mniejszy przedział ( a , 2 
bedziemy mieli wiec: 


dQ a, 


(7) dO 2 


+ b, . Sin © +a, . Cos 2 O b,.Sin3 O+ 
-+- a, . Cos 4 O+... 
wartości współczynników zaś: 


| Z 


TR 


2 
GË FSW WNE EE 
SE |i3-40» T ee OR 
R Lr 


2 2 


(8) z 
| p © Cos 20.40 
dG 


a to 


(8) że 


AIN 


| TO Cos 2 m O dO; bym=a, my1 =0 


2 


Wyznaczenie rozwinięcia (7) sprowadza się zatem do wy- 
znaczenia, drogą przybliżonych kwadratur, całek (8). 


W tym celu kreślimy najpierw krzywą, przedstawiającą 


w przedziale [5.5 funkcyę $S; trzeba bedzie wiec wyzna- 


czyć możliwie wielką liczbę punktów Ke Bereich czyli obliczyé 
dla danej szerokości e szereg wartości >> = < dla różnych © lub 8 


(Sin O= Sin e. Sin è). Obliczeń tych Geng Angot (właściwie dla 
d 3 
dal: 
stosując kwadratury przybliżone do całki (1) lub (2). Oto re- 
zultaty tych rachunków dla różnych szerokości geograficznych 
i różnych zboczeń słońca à, (p=0,7): (por. Tablicę na str. 50). 


ilości dziennej Q;, która tylko czynnikiem różni się od 


Jeśli liczby powyższe pomnożymy przez - , to otrzy- 


1 

1000 = 
i SKU 

mamy odpowiednie wartości dla „4, jeśli zaś pomnożymy 


2Cu 1 


Przez "oe": 1000 


to będziemy mieli ilości dzienne Q;. Mając 


e , wykreślamy odpowiednią krzywą i stąd 
szereg krzywych, których rzędne powstały od pomnożenia rzę- 


szereg wartości dla 


An 
dnych krzywej — odpowiednio przez SinQ, Cosic, it. d.; 


„Prace* Wydz. III. 08 15. Badania teoretyczne ilości ciepłą. 4 


50 — 


d 
Tablica Q; lub aa dla współczynnika transmisyi p = 0,7. 


ú Q : 
Table Qj ou 4. pour le coefficient de transmission p = 0,7. 


Szerokość geogr. ©. — Latitudes geogr. #. 
o° | 10” | 20° | 30° | 40° | 50° | 60° | 70” | 80° | 90° 


ND 
N © 
© 


556 
575 
590 
601 
607 
609 
607 
601 
590 
575 
556 
544 
536 


-+ 20° 


= — 
CO ND © 
© © © 


+ 40 


- declinaison du soleil. 


© 


Zboczenie słońca — à 
N 
© 
=] 


616 
621 
621 
617 
609 
596 
580 
559 
534 
506 
474 
457 
444 


653 
641 
626 
606 
582 
556 
526 
492 
455 
416 
375 
303 
357 


662 
634 
602 
568 
532 
491 
448 
405 
361 
314 
268 
245 
228 


682 
667 
646 
601 
554 
606 
456 
405 
354 
303 
254 
205 
159 
136 
120 


542 
506 
460 


514 
462 
393 


511 
454 
379 
237 
118 

34 
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zmierzywszy .pola krzywych otrzymanych, będziemy mogli obli- 
czyć wartości współczynników szeregu (7) ze wzorów (8). 

Dla Warszawy (w =5222) np. otrzymujemy następujące 
szeregi dla p=0,7 i p=0,8: 


p=0,8; Q;= 


(9) 


p=0,7;Qj= 


OC" " 


=i 


2Cu 


RO) 
i” 


(0,3842 + 0,3735 . Sin © —0,0261.C0s20+- 


+ 0,0045 . Sin 3 © + 0,0012 . Cos 4O +...) 


a? 
DH 


(0,3072+ 0,3105 . Sin©— 0,0289 . Cos2©+- 


+ 0,0046 . Sin 3 © + 0,1008 . Cos 4 ©+-...) 


eeh AE 
— (0,1223 + 0,1187 . Sin © — 0,0083.Cos20 + 
+0,0013 . Sin 3 © + 0,0004 . Sin 4© +...) 


Gm 0978--0,0988 . Sin ©— 0,0092. Cos2O+ 
+ 0,0015 . Sin 3 © +0,0008. Cos 4 O +...) 


Całkując (10), otrzymujemy wyrażenie dla ilości średniej Q: 


— . — C 
PE 0,8 , Q Wi c 
— 0,00415 . Sin 2 ©— 0,00048 . Cos 3 ©+ 0,00009 . Sin4Q©--...) 


(11) 


p — 0,7 ; Q == z 
— 0,00460 . Sin 2 © — 0,00047 . Cos3 © + 0,00007 . Sin 4 ©....). 


(0,11918 + 0,1223 . © — 0,1187 . Cos © + 


(0,09927 + 0,0978 . © — 0,09880 . Cos © — 


Jeśli oznaczymy przez te ilość ciepła, odpowiadającą 
okresowi, gdy słońce jest w półkuli północnej, Q—: gdy w po- 
łudniowej, przez Q; zaś roczną ilość ciepła, to, zakładając w sze- 
regach powyższych O= 7 i O=2r, otrzymamy dla p=0,8: 


p= 0,8; Qe = 0,6226 . - 
(12) Q— = 0,1458. 


Q: 


= 0,7684 .- 


106,840 kal. kg. 


25,020 kal. kg. 


I 


131,860 kal. kg. 


aja ola AD 


Dla współczynnika zaś p=0,7: 


C 
p =0,7; Q+ =0,5058. “= 86,800 kal. kg. 


(13) QE 


C 


.=0,1087.—= 18,660 kal. kg. 


Q:=0,6145. £= 105,460 kal. kg. 


gn va 


Porównywując wartość Q, dla p=0,8 z wartościami, otrzy- 
manemi dla Warszawy doświadczalnie, widzimy, iż wartość ta 
jest zbyt dużą; wartość Q;= 105,46 kal. kg. dla p=0,7 jest na- 
tomiast jedną z rzędu tych, którą zwykle otrzymuje się drogą 
doświadczalną; widzimy zatem iż średni współczynnik transmisyi 
dla Warszawy jest raczej bliższy 0,7 niż 0,8. Oczywiście, iż 
wartość 0,7 słuszną może być tylko dla większego przeciągu 
czasu, wogóle bowiem p ulega w ciągu roku znacznym waha- 
niom i dla mniejszych przeciągów czasu (np. dla miesiąca) re- 
zultaty otrzymane przy p=0,7 nie będą zadawalające; dla mie- 
sięcy zimowych np., jak później się przekonamy, należałoby dla 
p brać wartości bliższe 0,8. 

Jeszcze raz zaznaczamy, iż powyższe rachunki, jakkol- 
wiek mogą przy odpowiednio dobranym p, dać rezultaty zgodne 
z rzeczywistością, to jednak z powodu zbyt ogólnikowego tra- 
ktowania są tylko pierwszem przybliżeniem rachunków insola- 
cyjnych. 

Szeregi trygonometryczne, przedstawiające ilości dzienne 


AN 
Q; lub FE mogą służyć do badania przebiegu insolacyi dla 


danej szerokości +. Tak np. położenia słońca, dla których ilość 
dzienna Qj jest największą lub najmniejszą, otrzymuje się łatwo, 
przyrównywując do zera pochodną odpowiedniego szeregu dla Q;. 
Zauważmy przedewszystkiem, iż dla miejscowości strefy gorącej 
słońce jest dwa razy w zenicie w ciągu roku, należy zatem spo- 
dziewać się tutaj dwóch Maximów i dwóch Minimów dla natę- 
żenia średniego promieniowania lub dla Qj. Rzeczywiście, bio- 
rąc pochodną szeregu znalezionego dla równika, mamy następu- 
jący rozkład dla Maximów i Minimów: 


p= p = 0,8 p= 0,7 
GE A a EOT TEOSE e 


1 Max. 65°20' | 1010 68°10 737 69°10 617 
1 Min. | 170°10' | 888 16950" 629 169240" 517 


2 Max. | 273910" | 998 270720 | 726 269°10° | 607 
2 Min. | 3470 | 948 34780 | 672 341950 | 552 


— 53 


Dla miejscowości strefy umiarkowanej, będziemy mieli oczy- 
wiście tylko jedno Maximum, odpowiadające najwyższemu poło- 


żeniu słońca, t. j. dla © = A jedno Minimum dla najniższego 


2 


położenia słońca, t. j. dla © = —-r i to niezależnie od wartości 


współczynnika p (należy zwrócić uwagę, iż te położenia Max. 
i Min. nie będą odpowiadały rzeczywistości, gdyż są otrzymane 
w założeniu, iż p jest stałe w ciągu całego roku). 

Ciekawe jest wreszcie porównanie ilości ciepła, otrzyma- 
nych tego samego dnia pod różnemi szerokościami geogra- 
ficznemi. 

W dniu porównania wiosennego (o= o), największą ilość 
ciepła otrzyma równik; w miarę oddalania się od równika Q; bę- 
dzie stale malało; jest to zrozumiałe jeśli zauważymy, iż dla 
Ô—o długość dnia jest wtedy wszędzie jednakową, zatem na 
ilość Q; wpływa tylko wysokość słońca. W miarę jak 8 wzrasta 
i słońce oddala się od równika, maximum zmienia też swe po- 
łożenie, oddala się jednak od równika prędzej, aniżeli słońce; tłó- 
maczy się to wpływem na wartość Q; dwóch czynników, wyso- 
kości słońca i długości dnia; w miarę oddalania się od poło- 
żenia, gdzie słońce jest w zenicie, wysokość słońca zmniejsza się, 
lecz zato wzrasta długość dnia, ten ostatni czynnik początkowo 
przeważa i wpływa na znaczne przesunięcie położenia Maximum. 

Następująca tabliczka najlepiej przedstawi tę właściwość 
Max. Q; (zboczenie słońca równa się też szerokości, gdzie słoń- 
ce jest w zenicie, o odpowiednia szerokość, dla której zachodzi 
Maximum Ou: 

| Ge 8°; = 1280 
(14) = 16°; ę = 26°30 
| 0=e ; o = 43930 


e — otrzymano, kreśląc dla danych à krzywą przebiegu Q; dla 
różnych szerokości i szukając punktu najwyższego. 

Co się tyczy ilości rocznej Q+, jest ona największą dla 
równika i stale maleje w miarę oddalania się od równika ku 
biegunom. Nastepujaca tabliczka przedstawia wartości Q, (ilość 


PRASY DE 
średnią) dla różnych szerokości + i dla p = 1; 0,8; 0,7 (obliczenia 
według metod poprzednich): 


Tablica wartości rocznych ciepła, 
Table de valeurs annuelles de la chaleur. 


p p=l p=08 p=0,7 

0° 1,9182 1,3788 1,1458 
10° 1,8916 13540  1,1230 
20° 1,8132 1,2790 1,0544 
30° 1,6858 1,1624 0,9466 
409 1,5154 1,0078 0,8060 
50° 1,3116 0,8282 0,6488 
60, 1,0910 0,6474 0,4940 
70° 0,8855 ` 0,4827  0,3558 
80° 0,8002 0,3921 0,2690 
90° 0,7962 0,3738 0,2397 


Liczby powyższej tablicy należy pomnożyć przez 


(log É = 5,23453). 


8. Całkowita ilość ciepła otrzymana przez ziemię. 


W rachunkach obecnych przyjmiemy, iż ziemia jest ściśle 
kulistą. Łatwo zauważyć, iż wtedy obliczenie będzie można 
przeprowadzić tak, jakgdyby nie istniał ruch wirowy ziemi i jej 
ruch koło słońca, oba bowiem te ruchy nie wpływają wcale na 
kształt oświetlonej części, która będzie zawsze powierzchnią pół- 
kuli. Oznaczmy przez W całkowitą ilość ciepła, jaką odbierze 
ziemia; jeśli zastosujemy prawo Lambert'a, przyjmując pewien 
średni współczynnik p, to będzie: 

rr 
D EE E 

"8 

całka rozciągnięta jest na oświetloną część ziemi. Oznaczmy 
przez À, ę — współrzędne biegunowe, wyznaczające położenie 
punktu na kuli, będziemy mieli: 
(2) da = R? cosy dy. dA; sinh == cosy. cos À 
gdzie R oznacza promień ziemi. 


w 
I 


Zatem 


(3) W = BE 2 Kee 5 ET at 1 cos „cosA.dA. de" 
stad 
(4) W=W (0-0) 


całkowita ilość ciepła, jaką odbierze kula ziemska w danym prze- 
ciągu czasu, jest proporcyonalna do łuku jaki opisał promień wo- 
dzący ziemi w tym czasie. 

Zatem w ciągu roku ziemia odbierze ilość ciepła: 


de 8r CR? 


S fF (es .cosżę . COSÀ . dh. dọ 


Dzielac przez 4rR?, będziemy mieli średnią ilość ciepła, ja- 
ką otrzyma w ciągu roku cm? powierzchni ziemi: 


(5) We 


20 


S [= [> psec?.sech, cosżę . cos. dh. dọ; 


(6) Qu = 
D 6 

w przypadku p = l, t. j. gdy obliczamy dla powierzchni, leżącej 

na granicy atmosfery, będzie: 


(7) ee [> f Z cost. cosk.dh.d Ar., 
to jest: 
(7) Qu = 1,5708 . = = 269,5 kal. kg. 


Można obliczyć stąd, iż ilość ciepła otrzymana od słońca 

w ciągu roku zdolna jest stopić warstwę lodu grubości 33,7 me- 

tra, pokrywającą powierzchnię ziemi. Wartość (7) można otrzy- 

mać też łatwo bez żadnego całkowania. Zauważmy mianowicie, 

iż całkowitą ilość ciepła, wysłaną przez słońce ku ziemi, otrzy- 

mamy, mnożąc stałą słoneczną przez powierzchnię przekroju wal- 
ca, opisanego na kuli ziemskiej; będzie zatem: 

dW C 

Gore 


R. SEN 


dla średniej ilości Ou otrzymamy więc: 
dog CZE WIE gy WE 
dB ARR” BD = „4 MD (WST KARE 
wartość zgodną z (7). 
Przystąpmy teraz do obliczeń, uwzględniając stratę promie- 
niowania w atmosferze. Zadanie polega na znalezieniu wartości 
całki podwójnej (6), rozciągniętej na powierzchnię kwadratu, któ- 


rego bok wynosi Di można tego dokonać bądź przy pomocy przy- 
bliżonych kubatur, bądź też przy pomocy rozwinięć na szeregi. 
Zaczniemy więc od rozwinięcia na szereg, według wzorów 


rozdziałów poprzednich, funkcyi: 


wo 1 


9.eosh 
(OZ F (g) = | cosi. IT a 
J 
e loge E 
zaktadajac: Den Ee 
r Gs 
mamy: 
d F (4) = | cosk.e 7. di 
jak wiadomo: j 
5? Le "c Ce a, t a}. sin?p + a, sinte + a, sin*g +.... 


współczynniki a% Ay, a j t. d. są pewnemi funkcyami A, danemi 
przez wzory (17) rozdziału 5-go. Będzie zatem: 


[= [° 
(9) EC ar a,. cos . dh + 


o 0 


ksintp. | a,.cosk.d\+..... 


RE TT) 2. 


jeśli, podobnie jak w rozdziale 5-ym, założymy: 
T T 

1 ZDJ 1 

cosà 


( l 
om = | cosh. p Ou GH po dh; 
} 


D 
0 


fan 
In = E | secm-12.pseer, dh 
6 p I 
to otrzymamy: 
Ëss 
( H 2 | 
+= | aocosk.dh=L; n= | a,cosi.d=—gl,; 
J J 
1 
Y= | acosa. d =— gl Hol; 
J 
pi KC 
ATA a, cosà . dA = — 16 lı t 16 b — 4g h; 
(11) r 
Ze E EL 1 
Mira | a, cosà . dA = — 128 i + 128 2 64 ls t 3g4 143 
J 


© 


A ES BZ E 
5 — ] az cosà . dì = — ogg li + 956 l2 — 256 ls T 198 4 — 


D 
o 


1 
— 3849 ls: 


a BĘ) —= 


Obliczywszy zatem współczynniki je Tu "Ta... dojdziemy do 
rozwinięcia: 


(12)  F(5) = 1, .Ssin?ę +7. sintet ys sinp F... 


Wypada stąd: 


r 


i 
50 = + | (Usine) Frei. dẹ 


to jest: 


dW  4CR? (e E 
dose À [nt T) sine + (>t). sinte + 
J 


Dee E ] dę 


Posługując się wzorem: 


u 


k 
ge z 
Sinne. dọ = OEE 


otrzymamy rozwinięcie: 


dW EE | + 


dO ~ Zelt to) gy g GSA: Am: Ya) + 


35 , 63. Beat 
+98 (4 — YD * agg Us — 14) go fie 13) Te. = 
albo też: 


dW  2GzR*/l IZ Pal EC Ak 
do = c Mitt g" Tin qTo8!:' 256 4 t 


21 , 236 ea \ 
+1024 * 14336 67t) 


(13) 


EE e 
otrzymamy zatem dla całkowitej ilości rocznej, szereg: 


dar wwa SEE 
s r c 


1 1 5 

( to gii Figa" aah 
e? pa 21 yt \ 

* 256 = 10245 5) 


dzieląc zaś przez 4x R*, mamy ilość średnią dla 1 cm2: 


WA 1 l BIEŃ 
(15) On ne me en ta 


2214 Pire A 2 ` 
1024 ' 14336 55 ''''J 


dla p— 1 wiemy, iż I, = 1 zaś I, =l =... = Im = 0, wypada 
C 
G 


zatem y = 1, Y = Y, =... = 0, a więc Un = , jak być po- 


2 

winno. Zastosujemy teraz dane liczbowe; jeśli p — 0,8, mamy 

z rozdziału 5-go: I, = 0,7270; i, = 0,2204; I, = 0,0825 it. d.; 
stąd według (11): 


4 = 0,7270; w = — 0,1102; y, = — 0,0723; Y = — 0,0538; 
44 = — 0,0425; y, = — 0,0372; i t. d. 


Wstawiając te wartości w (15), otrzymamy ilość średnią 
ciepła, jaka dochodzi do ziemi w ciągu roku jeśli p = 0,8: 


DÉI 


(15) Qu = 1,0716 . — = 183 kal. kg. (p = 0,8); 
Jeśli przyjąć p = 0,7, KC wtedy: 


(17) Qu = 0,8662 . + = 148,6 kal. kg. (p = 0,7); 


DÉI 


Ponieważ p waha sie zwykle między 0,7 i 0,8 można więc 
przyjąć za średnią ilość ciepła, która dochodzi do om? w ciągu 
roku liczbę 165 kal. ke: porównywując z wartością (7) widzi- 
my, iż 104 kal. kg. czyli średnio 38% ciepła zostaje absorbowa- 
ne przez atmosferę. Jak widać z tabliczki wartości rocznych 
ciepła rozdziału poprzedniego, ilość średnia ciepła odpowiada 
mniej więcej szerokości 36°, 


E Ee 


9. Zastosowanie badań spektrobolometrycznych. 


W niniejszym rozdziale rozważać będziemy natężenie pro- 
mieniowania słonecznego oraz wpływ atmosfery, posługując się 
wynikami w dziedzinie badań nad rozkładem energii drgań 
w widmie słonecznem. Póki nie braliśmy pod rachubę wpływu 
atmosiery, rozważanie oddzielnych części widma było zbyteczne, 
prawo bowiem sinus'a i kwadratów odległości, stosowane w tym 
przypadku, słuszne było dla wszystkich rodzajów drgań. Z chwilą 
jednak, gdy w naszych rachunkach zaczynamy liczyć się z czyn- 
nikiem absorbującym energię drgań widma słonecznego, stoso- 
wanie ogólne praw Bouguera i Lamberta do całego wid- 
ma, przyjmując pewien średni współczynnik transmisyi, bez po- 
przedniego uzasadnienia, jest niedopuszczalne. Pomijać będziemy 
tę okoliczność, iż sama długość drogi atmosferycznej m, jest 
zależna od długości fali; trzeba jednak uwzględnić, iż atmosfera 
absorbuje w różnym stopniu różne rodzaje drgań, czyli, iż współ- 
czynnik transmisyi jest zależny od długości fali. Rozważmy bliżej 
ową stratę promieniowania w atmosferze. Stratę tę powoduje 
wiele czynników; po pierwsze istnieje tak zwana absorbcya 
właściwa, tutaj energia absorbowana zmienia swą postać, bądź 
przemieniając się w energię cieplną, bądź też w chemiczną, 
wielkość tego rodzaju absorbcyi zależy w wielkim stopniu od 
ilości pary wodnej w powietrzu. Drugą przyczyną straty energii 
jest to tak zwana dyfuzya czyli rozpraszanie się światła w atmo- 
sferze; ten rodzaj straty różni się od poprzedniego właściwością, 
iż energia stracona nie zmienia tu swej postaci; zresztą energia 
promieni rozproszonych nie jest całkowicie stracona, część ich 
bowiem dochodzi do ziemi, dzięki czemu widzimy przedmioty 
nieoświetlone bezpośrednio. Znaleziono, iż działanie chemiczne 
promieni rozproszonych jest równe 0,4 intensywności działania 
promieni bezpośrednich. 

Według powyższego, jeśli oznaczymy przez e», energię 
drgań X, a przez ei energię po przejściu atmosfery, to będzie, 
na zasadzie prawa Bouguera: 


s —(a--B)m 
(1) e= .E 


współczynnik a odpowiada absorbcyi właściwej, której główną przy- 
czyną jest para wodna, $ zaś jest współczynnikiem rozproszenia. 
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Według Ragleigh'a, ß jest odwrotnie proporcyonalne do Ae 
potęgi długości fali, wyrażenie przez niego otrzymane jest na- 
stępujące: 
2) BI) 1 
3N AN 

N — liczba cząsteczek powietrza w cm.*, p. — stała refrakcyi. 

Wzór ten jest zgodny z rzeczywistością przynajmniej o ty- 
le, (GB maleje, gdy zwiększa się À; powietrze rozprasza bar- 
dziej promienie o krótkiej długości fali, a więc niebiesko-fio- 
letowe, niż czerwone i infra-czerwone; tem zresztą tłómaczyć 
można błękit nieba. 

Opierając się na właściwości, iż promienie infra-czerwone 
w małym stopniu podlegają dyfuzyi, Radau (Actinometrie) po- 
daje wzór dla energii promieniowania na powierzchni ziemi; 
dzieli on mianowicie widma słoneczne na dwie części; dla jed- 
nej, zawierającej drgania infra-czerwone, przyjmuje współczynnik 


transmisyi p = 1 dla reszty zaś drgań p — Ss na zasadzie tego 


pisze Radau dla natężenia promieniowania Q- 


DIE 
Q = A, + A; (3) 

Dokładne poznanie praw absorbcyi i dyfuzyi energii sło- 
necznej, jest ważne w kwestyi znalezienia tak zwanej stałej sło- 
necznej, t. j. ilości ciepła, jaką odbiera cm? powierzchni na gra- 
nicy atmosfery; dokładna bowiem znajomość absorbcyi pozwoli- 
łaby ściśle z wartości obserwowanej dla natężenia promieniowa- 
nia znaleźć jego wartość na granicy atmosfery, to jest stałą 
słoneczną. 

Z przyrządów, służących do pomiarów natężenia promie- 
niowania, zasługuje na uwagę pyrheliometr systemu A ng stro m'a. 
Przyrząd ten posiada dwie cienkie blaszki platynowe, z których 
jedną wystawia się na działanie promieni słonecznych; powstała 
stąd różnica temperatur powoduje prąd, który kompensujemy, 
przepuszczając przez drugą blaszkę prąd z ogniwa, aż do wy- 
równania różnicy temperatur; z natężenia prądu kompensacyjne- 
go można sądzić o wartości natężenia promieniowania, 

Jak już wspominaliśmy, stała słoneczna jest to całkowita 
energia widma słonecznego na granicy atmosfery; jeśli więc 
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e, dh będzie energią drgań, zawartych między À i À + dh, to po- 
le krzywej, przedstawiającej rozkład energii w widmie ciągłym, 
czyli e, będzie właśnie ową stałą słoneczną; pole zaś krzywej, 
odpowiadającej przebiegowi e, t. j. energii widma na powierzch- 
ni ziemi będzie wyrażało energię promieniowania obserwowaną. 
Jeśli więc, mając krzywą energii obserwowanej wykreślimy (po- 
sługując się wartościami współczynników transmisyi, otrzymany- 
mi doświadczalnie dla różnych A) krzywą przebiegu energii dla 
widma na granicy atmosfery, to stosunek pól tych dwóch krzy- 
wych, będzie wyrażał stosunek między stałą słoneczną i natęże- 
niem obserwowanem skąd otrzymać można wartość tej stałej. 
Badania nad rozkładem energii w widmie słonecznem prowadzili 
ostatniemi czasy Abbot i Fowle w Washingtonie'). W ni- 
niejszej tabelce podajemy w skróceniu ich rezultaty, t. j. warto- 
ści e, (przypominamy, iż e, . dì będzie energią drgań À) na gra- 
nicy atmosfery i współczynniki transmisyi p dla powietrza suche- 
go (znaczy tylko wpływ dyfuzyi): 


02,60 0,80 140 | 16,3 106 | 21,0 | 210,5 
5,042 2,665 1,657 (0,89810,532 0,247/0,043 


| 
0,76 | 036 0,9 | 0,97 |0,97 0,9 | 0,9 


BR | 
À _|0u,30/07,35 04,40/0,45/05,50 
|| 
€ |0,539,2,684 1,338 6,057 6,062 
| | 
| px | 04 | 04 |0,54 0,64] 0,7 
I | | | 


Należy zauważyć, iż obchodzą nas tylko wartości propor- 
cyonalne do e;. 

Na fig. 2 przedstawiono rozkład energii widma na grani- 
cy atmosfery, jak również i widma na powierzchni ziemi (dla 
powietrza suchego) przy różnych wysokościach słońca. Wi- 
dać z rysunku, iż położenie Max. ea zmienia się zależnie od 
większej lub mniejszej absorbcyi powietrza. 

Położenie Max. ex, gra ważną rolę w teoryi promieniowa- 
nia ciał czarnych; według Wien'a, jeśli 7 — jest temperaturą 
bezwzględną ciała czarnego, zaś Am — długością wysyłanej fali 
o największej energii, to mamy: 


(3) Am « T = Const. = * 2930 


1) Abbot i Fowle: Annals of the Astrophysical Observatory of the 
Smithsonian Institution. Washington, 1907 i 1908. 


Prawo to, zastosowane do widma stonecznego, da pewne 
pojęcie o temperaturze słońca; w widmie obserwowawym na po- 
wierzchni ziemi Le leży zwykle między 01,5 i 0v,6, zależne od 
wysokości słońca; z rezultatów Abbot'a i Fowle'a sądzić nale- 
ży, iż dla widma na granicy atmosfery km = 0,47, co daje dla 
temperatury słońca wartość 7 = 6230%. Porównajmy tą wartość 
z wartością, którą otrzymamy przy pomocy prawa Stefana; 
prawo Stefana głosi, iż ilość ciepła, jaka w jednostce czasu 
wypromieniowuje przez 1 cm? powierzchni ciała czarnego o tem- 
peraturze bezwzględnej 7 jest proporcyonalną do 4-ej potęgi 
temperatury 7, zatem: 


(4) E=oT*; c= 768. 107% 

prawo to zostało potwierdzone teoretycznie przez Boltzman- 
ma. Jeśli 7 jest temperaturą słońca, to będzio („ — promień 
słońca): 


Ars 
a = stałej słonecznej = 2 kal. gr. min. 


„wypada stąd: 7 =5880? — tego samego rzędu co poprzednia. 
Naturalnie, iż znaleziona temperatura jest najmniejszą, jaką słoń- 
ce może posiadać, albowiem w rachunkach powyższych nie 
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uwzględniono absorbcyi promieniowania w atmosferze słońca; 
można zatem powiedzieć, iż temperatura powierzchni słońca jest 
większa od 6000". 

Wypadnie jeszcze zaznaczyć, iż dla ciała czarnego o tem- 
peraturze 7, Planck*) wyprowadził następujący wzór teoretycz- 
ny dla natężenia energii drgania A: 

(5) a=; C =5,29. 10°; C, = 14560; 
ell A 

Nie należy naturalnie się spodziewać, aby rozkład energii 
(5) zgodny był z rozkładem energii w widmie słonecznem, a to 
z powodu, iż warunki, pod którymi był wyprowadzony wzór (5), 
nie są spełnione w przypadku słońca. W każdym razie, trzeba 
zaznaczyć pewną zgodność (5) z doświadczeniem, wykazują to 
choćby badania Angstróm'a?). 

Rozkład energii (5) zgodny jest również prawem Wiena 
i Stefana. Otrzymujemy mianowicie, przyrównywując do zera: 

de. 

dh 
dowodzi to też, iż Max. e, jest odwrotnie proporcyonalne do 
5-ej potęgi długości fali jm. Aby z (5) wyprowadzić prawo St e- 
fana, zauważmy, iż energia całkowita widma będzie: 


-=0, iŻ An. T = Í (c; ©) = const.; 
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1) Planck: Vorlesungen über die Theorie der Wärmestrahlung. Leip- 
zig, 1906. 
2 K. Angstróm: Energie dans le spectre visible d’étalon Hefner. 
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zatem: Ed 
Ca I Log*z dz CG 
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o 
t. j. znane prawo Stefana. 

Niektórzy z autorów starają się wynaleźć wzory empirycz- 
ne dla wpływu dyfuzyi. Ciekawe są wyniki Angstrom'a w tym 
kierunku *): oparte one są na pojęciu widma fikcyjnego, którego 
wszystkie drgania mają jednakową energię; w tym celu zastę- 
puje się À przez ę(x), tak aby e,. dA = A . dx; następnie A ng- 
strom przyjmuje wzór empiryczny dla współczynnika transmi- 
syi: p = ax. W ten sposób dochodzi Angstrom do wyra- 
żenia wpływu dyfuzyi: 

KE pmó 
47-40: TE nma 
qə — stała słoneczna, p — średni współczynnik transmisyi, # i à— 
pewne stałe, wyznaczalne doświadczalnie: 


p= 0,93; n=0,26 
można więc napisać: 
ON 
1 +0,26. mó 


à — jest to pewna stała, która zależy od większej lub mniejszej 
intensywności dyfuzyi, można ją nazwać „gęstością dyfuzyi“. 
Dla energii W, absorbowanej przez parę wodną w powietrzu 
o ciśnieniu fr", Angstrom podaje następujący wzór, otrzy- 
many na podstawie wielkiej liczby doświadczeń: 

(7) W = gq,.0,1. (0,85) . (man 


stosując (6) i (7), będziemy mieli dla całkowitej energii, która 

dochodzi do ziemi, wyrażenie: 
(0,95)m% 

1 + 0.26. mò 


Wyrażenie to można stosować do wyznaczenia stałej sło- 
necznej ?). 


(6) i q — dg: 


Q = Go: — 0,1. (0,85)? . (Gmponnl ; 


1) Angstrom. Méthode nouvelle pour l'étude de la radiation solaire. 
Upsala. 


3 W. Gorczyński (l. ©). 


„Prace“ Wydz. HI. XN I5. Badania teoretyczne ilości ciepła. 


Przystąpmy teraz do obliczeń, będących celem niniejszego 
rozdziału, mianowicie do znalezienia rocznych ilości ciepła dla 
Warszawy. Oczywiście, iż w tym przypadku już nie będzie- 
my wprowadzali pojęcia ilości średniej ciepła, lecz użyjemy 
metody sumowań. Dla wpływu dyfuzyi będziemy stosowali war- 
tości podane przez Abbota i Fowle'a, zaś dla wpływu pa- 
ry wodnej weźmiemy wyrażenie (7). Naturalnie, iż wartości p 
dla powietrza suchego nie będą zupełnie słuszne dla Warszawy, 
otrzymane były bowiem dla Washingtonu. 

Na mocy powyższego trzeba będzie znależć szereg warto- 
ści dla natężenia promieniowania (nie uwzględniając wpływu pa- 
ry wodnej) dla różnych wysokości słońca. 

Wykreślamy w tym celu krzywą e przebiegu energii 
w widmie na granicy atmosfery i szereg krzywych p”e,, przed- 
stawiających przebieg energii widma po przejściu drogi atmosfe- 
rycznej m. Obliczamy następnie pola, odpowiadające tym krzy- 
wym w granicach od À (około 07,4) do À możliwie wielkiego; 
będziemy mieli zatem dla natężenia promieniowania d. wyraz 


następujący: 
KÉ 
| p” EA dA 


(8) äh an 


Otrzymamy w ten sposób szereg wartości żądanych: 


h = 90°; m= 1,00; q=0,818.q, 

h = 80°; m= 1,02; q= 0,810. o 

=70; m= 1,06; q= 0,800. qo 

h = 60% m= 1,15; q= 0,786. q 

(9) = 50°; m= 1,30; q= 0,7738. q 
h = 40°; m= 1,55; q= 0,760. q 

h = 30; m=2,00; q=0,714. a 

= 20°; m= 2,90; q= 0,635. qo 

h= 10°; m=5,6 ; q= 0,432. q 


EE EE 


Mając wartości (9) otrzymamy krzywą przebiegu natężenia 
d dla różnych wysokości słońca h. Posługując się tą krzywą, 
można będzie znaleźć dla q szereg trygonometryczny: 
(10) + = À, + A, cos 2n + A, cos 4h + À, cos 6h +. 


współczynniki otrzymamy, jak wiadomo, z kwadratur: 


2 q „gł WE 
Aua Se lg ANS s | q, DOS 2mh . dh; 


Mnożąc szereg (10) przez sinh, będziemy mieli wyrażenie 
dla składowych pionowych natężenia promieniowania: 


sinh.q 1 I A, Dez 
BA ENAR sinh + - 2 sin3h +.. 
qo qo 2 d t gd 4 
to jest: 


(11) i qu = 0,790. sin h — 0,035 . sin 3h — 0,008 . sin 5h — 
5 — 0,001. sin 7h—.…. 
Ciekawe będzie porównanie tego rozwinięcia z rozwinie- 


ciem, które otrzymamy dla wyrażenia Angstróm'a, kładąc 
è = 0,7; p = 0,93, mamy: 


ea E sinh = 0,797 . sinh — 0,047 . sin 3h — 
— 0,010 . sin5h —.. 


Stratę promieniowania wskutek pary wodnej, będziemy 
obliczać według wzoru (7), stosowalnego w granicach 5< f m<30. 
Postępując, jak poprzednio, otrzymamy rozwinięcie: 

Sin h . m0*% , 0,9" = 0,923 . Sin h + 0,027 . Sin 3 h -|- 0,005. Sin 5 h + 
+ 0,001 .Sin7h 

łącząc ten szereg z (11), będziemy mieli wyrażenie ogólne dla 

natężenia promieniowania słońca, gdy jego wysokość wynosi h 

i gdy ciśnienie pary wodnej jest f: 


(12) : Q= (0,790—0,092. #75) . Sin h — (0,035 +4- 0,0027 . f0-75) 


, Saski 
(0,008 + (0,0005 . {273} „5h — (0,001 +- 0,0001 . f0:75) „Sin 7 h 


SE 


W następującej tabliczce, będziemy mogli porównać natę- 
żenie promieniowania obliczone z (12) z wartością otrzymaną 
doświadczalnie: 

Q 
Obi.  Dosw. 
(kal. gr.) (kal. gr.) 
1912 r. AN ; 125.26p.; m= 1,46; fe 5,0%m, 0,87 0,82 
16/1 ; 12*.25p.; m=3,45; f= 1,7mm; 0,28 0,24 
AJNT ; 10*.48a.; m= 1,18; f=10,7"m; 1,04 0,87 
23j VIllz 7127: ; m=1,31; f= 9,6"; 0,88 0,76 
19/VI ; 10%. 4a.; m= 1,23; {= 7, mn: 0,98 0,96 
4IV ; 117. a.;m=l,5;i= 3,3mm; 0,82 0,89 


Możemy teraz przystąpić do obliczeń ilości dziennych cie- 
pła. W tym celu trzeba obliczyć pole, odpowiadające krzywej, 
której odcięte są kątami godzinnymi słońca, a rzędne natężenia- 
mi (9). Ponieważ mamy tabliczkę wartości natężenia g dla wy- 
sokości h, trzeba więc podać metodę wyznaczania dla każdego 
kąta godzinnego słońca odpowiednią jego wysokość. Bardzo 
prostą jest następująca metoda graficzna, podana przez Doleżala. 


Fig. 3. j 
Załóżmy we wzorze Sinh= Sing. Sin à + Cos 4 .Cos8. Cos T; 


Sin ę . Sin 8 = 


Sinh==y; Cost=x; V 
Cos » . Cos ê = m 


MG e 


otrzymamy wtedy linie prosta: 
y=mx+n 


której wykreślenie nie nastręcza trudności, jeśli zauważymy, iż: 
BB =mtn=_ Cos(4+8); (0B'=1) 
(CC)=m—n="Cos(ę—3); 


Dla każdego więc dnia, t. j. dla każdego 6, otrzymamy od- 
powiednią prostą, której część nad osią x-ów będzie odpowia- 
dać łukowi dziennemu słońca, zaś część poniżej osi x-ów łuko- 
wi nocnemu; punkt przecięcia się A odpowiada zachodowi 
i wschodowi słońca i daje jednocześnie wartość kąta godzinne- 
go przy zachodzie lub wschodzie słońca: Cost, =— (OA). Ła- 
two zrozumieć, iż dniu porównania wiosennego lub jesiennego 
odpowiada prosta, przechodząca przez O; gdy słońce jest w pół- 
kuli północnej, to punkt A leży na lewo od O, gdy zaś w po- 
łndniowej to na prawo. Otrzymawszy szereg takich prostych, 
nietrudno dla wszystkich dni odszukać wartości wysokości, a stąd 
i natężenia promieniowania, gdy zna się kąt godzinny; wystar- 
czy zresztą zauważyć iż < NOB'=t zaś < POB’ =h. 

Mając pole, otrzymane według wskazań powyższych, odej- 
mujemy od niego pole, przedstawiające energię absorbowaną 
przez parę wodną; dochodzimy w ten sposób do ilości dzien- 
nej ciepła poszukiwanej. Biorąc dla ciśnienia f wartości z roku 
1910, otrzymamy metodą sumowań następujące wartości insola- 
cyjne dla poszczególnych miesięcy (dla Warszawy). 


Styczeń 2,0 kal. kg. ; Lipiec 16,4 kal. kg. 
Luty 3,9 > 5 Sierpień 14,0 E 
Marzec 7,9 A d Wrzesień 9,8 Š 
Kwiecień 12,0 x Październik 5,8 A 
Maj 15,6 A 3 Listopad 2,5 G 
Czerwiec 17,1 ` ; Grudzień 1,4 A 


Wypada stąd ilość roczna Qr i odpowiednie ilości Q+- 
i Q—=: 
| Q:= 108,4 kal. kg. 1 cm? 
Q+:= 85,6 kal. kg. 
Q—:= 22,8 kal. kg. 


(13) 


ZW a 


Na zakończenie rozprawy niniejszej poddamy krytyce po- 
przednie metody traktowania absorbcyi powietrza. Przekonamy 
się teraz jak niesłuszne było tam przyjęcie pewnego stałego 
współczynnika p, będzie on bowiem w rzeczywistości pewną 
funkcyą wysokości słońca h, zależną też od wilgotności powie- 
trza. Jeśli q jest natężeniem promieniowania, to dla współczyn- 
nika pu, który będzie można nazwać średnim współczynnikiem 
transmisyi dla całego widma, otrzymamy: 


1 


m a \m 
gd, Di stąd: pu = E | 


posługując się powyższemi wartościami q, otrzymamy nastepu- 
pujące wartości dla p. przy różnej wilgotności f: 


Tablica wartości średnich współczynników transmisyi. 


h J= gomm f= 2mm f= gmm 
80% . . 0,818 0,710 0,656 
60° . . 0,811 0,711 0,657 
409 . . 0,839 0,754 0,709 


30° . . 0,845 0,773 0,741 
20? elle 0,796 0,762 
10° . . 0,860 0,811 0,790 


Widzimy, iż współczynniki wzrastają w miarę zbliżania się 
do horyzontu, znaczy to że dla wysokości mniejszych strata energii 
na jednostkę długości drogi m jest mniejszą niż dla wysokości 
większych. Łatwo będzie na zasadzie tego wytłumaczyć pewne 
niezgodności rezultatów rozdziałów poprzednich z doświadcze- 
niem. Jakkolwiek przyjęcie p=0,7 doprowadziło w rozdziale 
7-ym dla ilości rocznej Q, do wyników zadawalających, to je- 
dnak dla okresów czasu mniejszych istnieć będzie pewna nie- 
zgodność z rzeczywistością; tyczy się to mianowicie stosunków 
mamy bowiem: 


x 
(a; 


p=l; Al 299 


= 


Lee Ee 


y Qe 105,8 


tte 50 


= 4,27 


EE 
p=0,; GER 18,7 4,65 


według zaś metod spektrobolometrycznych wypada: 


s ren GENRES 


GËT 22,8 


liczba ta jest najbardziej zbliżoną do liczby otrzymanej ze, śre- 
dnich pomiarów w r. 1905—1910: 


Q+e — 3,60 


m 


ke —£ 


Główną przyczyną tej niezgodności jest to, iż przyjęto dla 
współczynnika p tę samą wartość w zimie co i w lecie, wskutek 
Q 


+e 5 > ń 
— za duże w porównaniu 


z doświadczeniem. W zimie bowiem, mamy dwie przyczyny, 
które wpływają na powiększenie współczynnika p; po pierwsze 
absorbcya pary wodnej jest mniejszą (średnio wtedy f= 2mm) 
powtóre zaś, wysokość słońca jest niewielka, co, jak widzimy, 
z tabliczki średnich współczynników transmisyi, również wpływa 
na powiększenie współczynnika p. Ostatecznie zatem, chcąc pro- 
wadzić rachunki insolacyjne według metod rozdziałów poprze- 
dnich należałoby dla miesięcy zimowych przyjmować p raczej 
bliższe 0,8 dla miesięcy letnich zaś brać p mniej więcej w gra- 
nicach 0,66 i 0,75. Wogóle zaś wypadnie zaznaczyć, że ścisłą 
metodą może być tylko stosowanie badań spektrobolometrycznych. 


czego Q—: było za małe, a zatem 
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Recherches théoriques sur les quantités de 
chaleur, reçues par la terre, en tenant compte 
des pertes du rayonnement dans l'atmosphère. 


l. Considórations générales. 


En premier lieu de nos recherches théoriques nous allons 
faire le calcul en supposant que toute la quantité d'énergie que 
le soleil nous envoit arrive à la surface de la terre, autrement 
dit nous ferons le calcul pour la surface horizontale à la limite 
de l’atmosphère. Après cela nous allons modifier le calcul d'in- 
solation, en tenant compte de l’absorption dans l'atmosphère et en 
prenant comme base la plus simple expression pour les pertes 
du rayonnement, notamment la loi de Lambert. A la fin de nos 
études nous nous basons sur les résultats des recherches dans la 
spectrobolométrie. 


En désignant par A la hauteur du soleil, » — distance du soleil 
à la terre, nous avons pour la quantité dq de la chaleur, reçue 
dans le temps dé, l'expression: 


k 
(1) dq =% Sinh.dt 


C čest une constante telle que > exprime la quantitć de la 
chaleur reçue dans l'unité de temps à la limite de l’atmosphère; 
z réduit à la distance moyenne s'appelle „constante solaire“, sa 
valeur est, d’après les mesures de M. G or cz y ń ski, égale à 2,05 gr. 
cal. min. I cm2. La quantité diurne de la chaleur on obtiendra 
en intégrant (1) dans les limites relatives au lever et coucher du 
soleil. Cette quantité, que nous calculerons en considérant les coor- 
données du soleil constantes pendant la journée, sera une fonc- 


Kë E 


tion évidemment discontinue de ces coordonnées car la longitude 
du soleil ou son déclinaison changent brusquement de jour en 
jour. Pour calculer la chaleur reçue dans l'intervalle plus grand 
de temps, correspondant aux longitudes du soleil (©,, ©), il n'ya 
qu’à faire la somme de quantités de la chaleur pour les journées 
particulieres, de telle sorte qu’on peut écrire: 


To 
©: ©: P: 
2 Ez 
(2) Ga Ss (HE 


f (7) — intensité du rayonnement, % la moitié de Parc diurne du 
soleil. Nous allons dans les chapitres suivants simplifier le 
calcul (2) à l’aide de la conception de la valeur moyenne de la 
chaleur. 


Il. Quantité de chaleur pour l’équateur et les pôles. 


On a pour la quantité diurne Q; l’expression: 


(1) O (EZE) 
0 
si e exprime la latitude géogr. du lieu, à — déclinaison du soleil, 
T— langle horaire, on a d’après légalité connue: 
(2) Sinh = Sin #.Sin à + Cos 9. Cos à. Cost 


en remarquant aussi que: 


86164 d 4 e 
t— DE, OU Hrs: op Puisque la journée stellaire 
a 86164 sec. moyennes, on aura donc: 


(3) Q = ŚCIE (Sin e. Sin ò. ze + Cos ẹ . Cos è. Sin T,) 


Tọ nous calculerons en ćgalant dans (2) k = 0; il en sera: 
(4) Cos T, = — tg ọ .tgò. 
Pour l'équateur 9 =0, ger par conséquent, la quantité 


diurne sera: 


Di 
(5) Qe = SSC Cos à. 
et pour les pôles: © = zo = T, donc: 


(6) Qp = SC e. Sinè. 


En se servant de la formule (3) on a aussi construit la ta- 
ble des valeurs de Q; pour les différentes latitudes et les épo- 
ques de l’année (chapitre 2-me du texte polonais). Pour avoir 
les quantités Q; en gr. calories il n’y a qu’à multiplier les nom- 


bres de cette table par tel que: 


EL 
1000 ro 
a 
log —— = 2,96105 
cette valeur est donnée pour valeur moyenne der si l’on désire 
l'exactitude plus grande il faut appliquer la formule: 


(z) ES - e . cos (© —w)]? 


r Fare? 
e — 0,0168; „ — longitude du perihélium, ou bien: 


D 12 
7) L — 1,00014 + 0,03660 . Cos $ + 0,00071 . Cos 2$ 1... 


$ — l'anomalie moyenne. 
Calculons maintenant les quantités de chaleur, reçues dans les 
intervalles plus grandes de l’année sur le pôle. On a: 
dQ Gars 
Sin à 
dt > Sin 
mais en remarquant qu'entre la longitude © et la déclinaison à 
existe une relation sin = sin e. sin © (e— inclinaison de l’éclip- 
tique, égale à 232730”) on a: 
dQ CE 
Z À Sine. Sin © 


D’après l'égalité connue de cinématique: 
p g q 


EE 


dO 
r? dt =€ (constante des aires) 


(ra 
donc: E — Sine . Sin O; 


d© c 
par conséquent l’expression voulue pour le pôle sera: 
(8) Q® Ee TEE Ee EN 
Oo C 


en prenant pour l’origine I'ćquinoxe ©,— 0, on a alors: 
(9) Qa = 2 Sine . Sin? > 


De lå résulte la quantité annuelle Q, reçue à la limite de Pat- 
mosphère en substituant © = 27: 


Qr = 2 A Sin e = 0,79620. CS 


ou bien en calories: 
Q: = 136,626 kal. kg. 


en remarquant que: 
log g = 5,23453 


Considérons maintenant le cas de l'équateur. Au lieu de 
faire les sommes de la chaleur indiquées par (2) du chap. 1 nous 
allions calculer la chaleur moyenne, indiquée dans la suite. Con- 
sidérons pour cela la quantité de chaleur reçue par une bande de 1 
cm. de largeur située tout le long du périmètre de Péquateur. 
Désignons cette quantité par W; ce sera la chaleur, reçue à Pin- 
termédiaire de la moitié éclairée de l’équateur; il est évident que 
W ne dependra que de la position du soleil et sera une fonction 
continue de leurs coordonnées. Si l'on divisera la quantité W par 
le périmètre on aura par conséquent la quantité moyenne de la cha- 
leur reçue par l'unité de surface de l'équateur dans l’époque 
correspondante. Nous allons calculer W; on a, en désignant 
par X, l'arc compris entre le point considéré et le point où le 
soleil passe par le méridien 


ként y EE 


(10) | Sinh = cos. cos À 
|  da=Rdh; R—rayon de la terre 
on a donc: 
ZĘ J sin TEE $ ee) Béi 
Tee das OSÓ.COSA. 
eq bd 
T2 
ć dW _ 2CR À 
de là: TIG. = pa 1 cos 8; 
en dćsignant par Q la quantitć moyenne: 
1 
WAĆ 27R W 
AQ 5 
Il en sera: EE cos 8 (11) 
en éliminant £ d’après légalité r? ER =c On a: 
ACY 1 C 
7X = — —. 8 
(12) dO z c CoS 


mais: Sin à —sine. sin—©; cos 8 = V 1 — Sine. Sin? O; 


on aura finalement la quantité cherchée: 


„O 
KG 
z e | V1-Sine.Sin© .d© 


Oo 


(13) Q= 


T ; 3 
en posant © =0; O=—, nous arrivons à la quantité annuelle, 


reçue à l’équateur à la limite de l'atmosphère: 


(14) Ch) EE do 
0 
C.— à. — d. intégrale elliptique complete de 2-me espèce, que 
lon désigne habituellement par E (e). 
En se servant de la table des valeurs de cette intégrale 
(„Tablice Matem.-Fiz.“ du Prof. Witkowski) nous avons: 


E EE 


4 C ; C 
(15) Qr EE E (23927') = 1,9182 z 
ou en calories: Q, = 329,170 kal. kg. / cm.?; cette valeur est bien 
acceptable si nous la comparons avec la valeur, obtenue par la 
méthode des sommes: 

Qr = 332,640 kal. gr. 
cela confirme nos méthodes. 

On pourra aussi représenter la quantité (13) par une série 

trigonométrique. En effet d’après la formule de Newton: 


16 SE sin? sin? © — À sin‘ e. sin* © — ) 
( OTS c 9 e. 8 ne. OCET 
il n’y à qu'à transformer cette série. Considérons plus généra- 
lement: 
(17) S =a, + a, sin? x + a, sint x + a, sin x +.... 
on a: 

Sin 2x= L l cos2x 

3 

i E 1 cos2x+--cos4x 
(18). GE 8 

Site Dee 2 +Š Os 4x — | 056 

in *x = 7g COS2X FE C 30 x 


en substituant ces valeurs dans (17), on aura: 


(19) S = A, + A, cos 2x + A, cos 4x +... 
ou: l 3 5 
A. ao Ty 4 Fog 42 | 16% 
fi E AT t 
(20) 1 9 "i 9 ie 39 43 
] 3 
AA g 42 116431 


(regardez la formule 33 du chap. 2-me du texte polonais). 


EE reż 


A l’aide de (20) nous aurons pour (16) le développement de 
Fourier: 


dQ l 
SK SC E + À, cos 20 + A, cos 4O +...) 
EPEN E NS. PR 
EA g Sin? e Së ECH 
(21) A= 4 Sinte+ emgeet, 
A. = TERS ee 
2 64 e 


en appliquant les données numériques: 


O = (0,30530 + 0,01318 Cos 20 — 0,00014 . Cos 40 +.) 


et en intćgrant on aura la sćrie voulue: 


(22) 


(23) Q= € (0,30530 . © + 0,00656 . Sin 20 — 0,00003 . Sin 4© +...) 


en posant © = 27, nous aurons une UE concordante avec (15) 


Qr = 1,91824. = 


Ill. Recherches des quantités de la chaleur pour une latitude 
quelconque y. 


Nous allons commencer tout de suite le calcul de la quan- 
tité moyenne de la chaleur dont la définition du paragraphe 
précédent s'applique aussi dans le cas actuel. Nous considé- 
rons donc une bande de 1 cm. de largeur, étendue tout le long 
d’une parallèle. Si À a même signification que précédemment, 
nous aurons de même: 

f Sin À = Sin y. Sin à 4- cos 4. cos a. cos À. 


(1) 
|  ds=R. cosg. da. 


ES 


En prenant l'intégrale étendue le long de la partie éclairée 
de la parallèle: 


À 
= | f (Sin o Sind+cos ẹ . cos ò . cos À) R . cosy. dà. 
‘0 
où cos ły =— tg ọ. tg è 


À — c’est d’ailleurs l’angle horaire. On a, en intégrant: 


De = PCR 2081 d (Sin e, Sin 8. Ag + cos p. cos à. Sin Aal 


pour la quantité moyenne Q nous avons de même: 
E Ae 
RS 2m R.cosyę' 
puisque 2m R.cosę est le périmètre de la parallèle. 


Nous arrivons donc A la formule: 


AN 


(2) REZ (ny. Sin 0.4, + Cosę. Cos8. Sind) 


ul 


en remarquant que Sin —Sine.Sin©, nous aurons finalement: 


(3) AA (E Sine.SinQ.4,+- 


SS S Sin? 
y 1 cosie e . Sin? ©). 


Cette expression laisse aussi se développer en série de F o u- 
rier, Il faut se servir du développement: 


4 = + are Sin (tg ę .tg 0). = 


T || 
= + (87.180) + -g (tge. te) + 
nous avons ainsi: 


,=-5 + Sine. SinO tg g + Sin? e. Sin Greet tg al 


l- Sin° s . Sin* © TE tg? pHi tg > +: 


RORY. R 


de la nous obtiendrons: 
tgv. Sin e.SinQO ds a Sin © + a; . Sin? O + a; . SintO p... 


a, = 5 «tg. Sine 
a = tg?e. Sin?e 


ER LE o tę + - ee l Sint e 


(formules 10 du Chap. II du texte polonais). 


Cos? 
me dans le chapitre 2 (formule 13 Chap. II du texte pol.) De 
cette façon nous aurons le développement: 


(4) CN. (M, +M,. Sin © + M, Cos 2 O + M; . Cos 4© +.) 


et en intégrant: 


Le développement de YA IK Sinte Sin? © sobtiendra com- 


(5) = Č M, +M,0 — M, Cos O + - E Sin? O+ 


EC 4 Sin4 O+...) 


Désignant par Q+: la chaleur correspondante à la position 
du soleil dans l’hémisphère Nord et Q—. dans l'hémisphère 
Sud, on a: 


Q+e = S KEEN 


(6) x 
Q—: = ve (xM,—-2M,) 
et la quantité annuelle: 
Qr = M . 2x M,; 


c 


„Prace* Wydz. III. M 15. Badania teoretyczne ilości ciepła. 6 


ou 
C 1 3 15 
T Q= Ek Los EN + gas + ju te) 
i 1 ëmze 3 Sinte 
NE Cosżę Dose 
De (6) résulte la relation intéressante: 
Ć 
Que 7 Q—e = "e 4M, 
c.— à. — d. 
C 
(8) Q+:—Q—=2 —. Sin. Sine. 
En appliquant ces résultats pour Varsovie (p=52,22) (re- 


gardez les formules 21 et 22 Chap. 4 du texte pol.) nous aurons: 


(9) = © (0,2007 + 0,157. Sin © — 0,0066 . Cos 20 +...) 


ép 


(10) Q=— (0,1573 -+ 0,2007 . ©— 0,1573. Cos © — 


— 0.0033 . Sin 2© +...) 
de la: 


Q-:=0,9451. = 162,2 cal. kg. 


eng 


Ee Ee EE 


| 1= 1,2610. € = 2164 1 
du 


Il y a une relation simple entre la dérivée iG et la quan- 
tité diurne Qj de la chaleur. 

En effet dans Pune et dans Pautre nous avons les mêmes 
intégrations à effectuer: 


al (Sin 4 . Sin à + Cos ẹ . Cos à. Cosi), dà 


d 


To 
Qj Sy f (Sin . Sin 6 + Cos. Cos è. Cosch de 


u 


a RYSA 


il résulte de là la relation fondamentale: 


l dQ  2rCu dQ 
(12) Qj = Ham. T 


il faut remarquer seulement que pour = © varie ďune façon 


continue ce qui n’a pas lieu dans le cas de Q;. On déduira 
donc facilement de la série (9) Pexpression suivante de la cha- 
leur diurne pour Varsovie: 


G 
(13)  Q;—0,001 SC (3,443 + 2,698 . Sin © — 0,113. Cos2 © + 


+ 0,007 . Cos 4 O --...); 
Tous les calculs précédents sont applicables pour o < 90° 
— e= 660838; pour e plus grande que cette valeur, c. à. d. pour 
la zone polaire, il faudrait faire le calcul sćparement pour les dif- 
férentes époques de l’année. 


IV. Les pertes du rayonnement dans l'atmosphère. 


Pour tenir compte de lPabsorption dans l’atmosphère nous 
allons prendre comme loi exprimant la perte de l’enérgie, la loi 


de Lambert: 
1 


Sin h 
q go -P 
ge — constante solaire, h — hauteur du soleil. 

La loi de Lambert n’est quapproximative pour I atmos- 
phère; nous allons l’établir. On a, en effet, d’après la loi de B o u- 
guer pour l'intensité q du rayonnement, après la traversée d’une 
couche o homogène: 


—ko 
(1) q=q,:€ 


k — constante spécifique. Puisque l'atmosphère n’est pas ho- 
mogène on ne pourra appliquer (1) que pour une couche in- 
finiment mince; on a alors pour la perte — dq de la chaleur dans 
cette couche l’expression: 

(2) da=—K(H).q.ds 


ELE 
K sera maintenant une fonction de l'altitude H. 
En supposant que le chemin parcouru par les rayons lumi- 


neux est rectiligne on a: c= -zi> par conséquent: 
JJ 


111 11 


dg a 
de là: A 
re EEN 
Log, (a) Sint | K(H). dH 


H est aussi grand que possible; q — rayonnement à la limite de 


l'atmosphère égal à 2 cal. gr. min. En posant: 
H 
fx. au 

(3) DRE 

nous aurons la loi de Lambert: 


1 
Sin h 


(4) q=qo-P 
h—s'appelle le „coefficient de transmission“. 

Nous allons maintenant généraliser la formule (4), en suppo- 
sant que le chemin des rayons lumineux est curviligne et en te- 
nant compte de la sphéricité de la terre; en admettant que K=ap 
, — densité de la couche à la hauteur H on a: 


H 
al lf Ze D 
Log, Ex | dé o da 
o 
posons: 
H H 
BR EI IO 
©) PoS WA mA 
IER 
0 
il en sera: 


(6) q=4.P ; (P&1; m2!) 


m qui pour h=o se réduit à 1 s'appelle „la masse atmosphéri- 
que traversée“; m, comme l’on verra, diffère d’ailleurs peu de 
1 
Sin h ` 
(R +H) w. Sin i = Const. ź langle d'incidence pour la couche con- 

siderée, u — indice de réfraction. On aura: 


Il n’y a qu’à calculer m; d’après la formule de réfraction 


H 
b CR, 
o m=fogu | 1/1 (Bis. Vatnss 
e dH | Le (R Sn Sin? z 
H o x 


en supposant que le chemin est rectiligne (# =), on a pour m 
la série suivante: 


H 
fond 
m =secz— secz.tg*z żę RA 
R IER 

(0) 
si Ton admet que p diminue dans une progression géométrique, 
il sera: 
(8) m= sec zk. secz. tg?z + (l, — RÓB lo tg" z 

2R 38 2R? Cosi: 


lọ — hauteur de l'atmosphère, supposée homogène et de densité po. 
On voit que l'expression (8) de m pour la distance zénitale pas 
1 


grande diffère très peu de => 


Sinh? dans le calcul suivant nous 


prenons toujours RE Il y a aussi à noter la formule de 


Laplace: 
Refr. 


— K. Cosh 


que l’on pourra établir facilement en prenant les formules de la 


m 


réfraction astronomique. 


GER dE 


V. Calcul de la quantité de chaleur pour l’équateur en tenant 
compte de Fabsorption atmosphérique. 


En nous basant sur la loi de Lambert nous aurons: 
a 
1 


Sin h 
HS Sin h p 


en remarquant que pour l’équateur sin h = Cos ò.Cosqt, on a la 
quantité diurne de la chaleur: 


D 1 
a Cos à. Cos T 


AO 
(1) Q= Se | Cosż. ge „dz 


“o 
et pour la dérivée de la quantité moyenne Q: 


WÉI 1 
gek Cos s . Cos À 
(2) a E ~ ] Cos 8. CosA.p . dà 
j 


De même comme Mr. A ngot nous allons chercher le dévelop- 
dO "KŚ ` A 
pement de Qj ou gg en série de Fourier, seulement au lieu 


d'utiliser la méthode d'interpolation d'A ngot nous allons repren- 
dre la question théoriquement. L'expression (2) pourra s'écrire 


aussi: Sec À 


= / Was 
x |  !-Sin*eSin*O.p 
SL 


D 
O 


„Cosa. dÄ 


Toute la difficulté du problème consiste dans le développe- 
Sec À 
ment de la fonction V 1 ginge. Sin ©» OÙ; en posant x= 
P 


Sin e. SinO; a= Sech. log, +: 


A RETIRE 


(4) Ve ; a>0, 


c’est une fonction holomorphe à l’intérieur du cercle de rayon 1; 
nous démontrerons que la fonction (4) sera représentée par une 
série de Taylor: 


Re e 
dure 


(5) e =a Bay XF dą X! Fay x +... 


on a en effet, d’après la théorie générale de fonctions analytiques, 
l’expression suivante de la rente: 

xn+t p f(z).dz 
(6) ee 
C Ki 


Vi? 


x— point intérieur du cercle dans le cas actuel f(z)=e 


Le module de l'intégrale le long d’une courbe est moindre que 
le produit de maximum de la fonction par la longueur de cette 
courbe, par conséquent: 
R |< A R RE Dëse (z)=1 
* "IR Maxje Ti Zoe 
On trouvera facilement que pour le cercle de convergence: 
(z=ģtin: 


DES e EE E 
Maxle His Max e ? Lt Vi 


NK St? 


on a donc la limite cherchée: 


j v 24+2ÿ 2 
|yln+1 


: e 
SIZE 


(7) . IR, 


puisque |x| <1, Ral tend vers zéro lorsque n tend vers l'infini, 


donc la série (5) est convergente, ce qui fallait démontrer. Il faut 
maintenant reprendre le calcul des coefficients ag a, a, a}... 

Nous ferons cela de la manière suivante; considérons pour 
la simplicité le développement: 


a 


ae 
vea A E Ee 


a 


À dy re ki 1 À -8 
de là: LEI HE — Z — — 0 1—x £ 
De 2 (1 — x} 9 y ( ) 
d Lane ARS E EE ere 
g otactae+...||(1+getgeryge + Tnt 


en multipliant les 2 dernières séries et comparant avec la série: 


d 
Y = ay F 243X T JA ES ch 
on aura: dx 
R 1 Ee 
za a. purs 
4, —€ a De Ge dz — 4 lo % raj] 


(regardez la formule 16 du Chap. 5 du texte pol.). Les formules (8) 
permettent de calculer les coefficients, on déduit de là: 


e > — 1 
4, = € 


(9) 


ae~? 
a =— 6 — 2) 


formules 17 du texte pol.) 
En appliquant le dćveloppement considerć, on a: 


zc d Sat, P 
C dO —VI=-sinte. sim O. | | a,. COS À. dÀ + 
O 


sin? e . sin ©. | i a, Cos À. dA + sin‘e.sin‘©. | “ap. cos h.dX-|... 


O 


SE 


nous arrivons finalement au développement: 


(10) D — to + ni Sin © +. int © +. sin O+ SÉ: 
Ou on a pose: 
To — lo 
SET 
sin*e 


(1) 
| 2 ~-z ++) 


(formules 19 du texte pol.) 
les quantités I sont données par les intégrales: 


(12) L=| "psk „cos. dh; In= log? 1 p secta pee dis 
o Pa 
O 
de (10) on déduit facilement le développement de Fourier: 


(13) TO = (A, +A, . Cos20+ As. Cos40+....) 


les coefficients sont les suivants: 


1 3 E 
Art ht g fat ZYWO dus 


16 


do PAIR TMZC E 
(14) dl 9 fi 9 la 39 3 16 4 
Lt T 


À, = 


Quant au développement de Qj on a d’après la formule 
_ ZC dQ, 
de 


un Q =Ç" (A, +A, .Cos20+A,. Cos4O +...) 


Get ul 
en intégrant E nous arriverons à la série cherchée: 


(16) Q == el. SÉ !$in20+ $ ZAŚ. | 


qui représente la quantité moyenne de la chaleur reçue dans 
l'intervalle (O, ©); pour la quantité annuelle de la chaleur on 
aura l'expression en posant © =25r dans (16), donc: 


90 I 1 a E OE 
(17) Gage KO | 
Comme exemple numórique prenons p=0,8, alors: 
l, = 0,7270; I, = 0,2204; I, = 0,0825 
de là: Yo = 0,7270; yı = — 0,0780; ya = — 0,0081; 
alors nous avons 
dQ 1C 
EE .6894 + 0,0380 . Cos 20 — 0.0004.Cos 40 +...) 
Q= 
On a aussi: 


= (0,21944. © + 0,00604 . Sin 2 © — 0,000083 . Sin 4 ©—....) 


al 


57 (0,6894 +- 0,0880 . Cos 2 © — 0,0004 . Cos 4 © +...) 


de là on aura la quantité annuelle: 
Q, = 1,3788. = 236,610 cal. kg. cm? 


Mr. Angot a obtenu pour la méthode des sommes 236,1 cal. kg. 


VI. Calcul de la quantité de chaleur pour le pôle. 


Ici on n’a pas besoin de la conception de la chaleur moyenne, 
puisque le mouvement diurne n’influe pas sur l’insolation. Pour la 


quantité diurne on a tout de suite l'expression: 
1 


(1) Q= "SA Sinz. ont, 5— inclinaison du soleil. 


et pour la CAE reçue dans l'intervalle correspondante (0,0): 


Yo (fre 
C 9 Cosc e. Cosc © 
(2) Q=< [Sine.Sin©.p 
0 
il ny a qu'à développer cette quantité par les formules du 
chapitre précédent; on a: 


a 


Cosc :.Cosc() CS Vi gese 
p =e = ao + a, Cos? © +a, Cost O +.... 
1 l 
en posant a = Sin: loge p' 


Il faut remarquer qu’on dćduira tout de suite la quantité an- 
nuelle: 


C p2, 
QEA | Sine. (a+ a, Cos? © +...) . d (Cos ©) 


0 
donc: 


ZE. m DEE 
(3) Qr Se WE 


pour avoir le développement de Fourier, on trouvera: 


Cosc e.Cosc() 
D = Ba + f, . Cos 2 O +f,. Cos 40+ $,. Cos6Q0--... 
par suite: 
dQ_C i SAN i H -P D 
(4) jj =< Sine (fo 3 |Sino+ me SCO 


et en intégrant: 
(5) Q= © Sine (AA. Cos O+ A, Cos3 O+ As. Cos5 O+...) 


A, A,... sont donnés par les formules (9) du texte polonais. On 
appliquera les résultats pour p= 0,8; les valeurs de a, a, a, sont 
données par la formule (10) du texte polonais; nous aurons: 


= > (0,18691 — 0,20117 . Cos © + 0,01137 . Cos 3 © + 0,00221 . 
Cos 5 © + 0,00039. Cos7 O+....) 
d’où résulte la quantité annuelle: Q, = ©. 0,37382 — 64,15 cal. kg. 


A SODA 


Vii. Calcul des sommes de la chaleur pour une latitude quel- 
conque % en tenant compte des pertes de la chaleur. 


En appliquant les formules précédentes, nous aurons lex- 
pression "ien 


(D Q= RE * fine. Sin è + Cos4 . Cos à. Cost). 
0 i 
Sin ę. Sin à + Cos. Cos à . Cos t 
dr 


et pour la quantité moyenne Q: 
LR OMS 
(2) Ee f Sing . Sin à + Cos ẹ . Cos à. Cos À). 


0 o 
1 


Sin o. Sin à + Cosy. Cost. Cosi 
. P „dh 


De même que précédemment nous allons chercher les dé- 
veloppements de Fourier pour (1) et (2). Mais la méthode 
précédente, qui s’appuye sur le développement de Taylor, pré- 
sente dans le cas actuel une difficulté considérable; nous utilisons 
une autre méthode en nous appuyant sur la théorie des séries 
trigonométriques. En effet chaque fonction satisfaisante aux 
conditions de Dirichlet, peut être développée dans l'intervalle 
(a, b) en série: 


F (x)= (t)= | (a, Cos t+ b, Sin t)+ (a, Cos2t+-b, Sin2ty+... 


EE SAM 
e Ire AB) = few, dt; 
am = — [Cosmt. om. Ar be (ënn A 


Leg =t 


d 


Es l 
Appliquons (3) à la fonction périodique ; "w 


SE EE 


Si l’on remarque que Melon dO }: 
on aura: 


(4) FA = %+b,.SinO+a.Cos2©+b,.Sin©+ a, . Cos4O+... 


27 
aaea TR - Cos2©. dO; byn= amer 0; 


Al 
Il faut donc obtenir une série des valeurs de Pe à l’aide 


de la quadrature approchée; de là on construira les courbes Es 
) 
R. Sin © etc. et l’on trouvera graphiquement les coefficients 


10 
a, b, & etc. (regardez la table de valeurs de E? texte polonais). 


De telle façon on trouvera les développements voulus de Q;, 
dQ 
d© 
7 du texte polonais (formules 9, 10, 11, 12, 13). On a conclu 
que la valeur p=8,7 donne en moyenne pour Varsovie des résul- 
tats du même ordre que les valeurs observées. Il y a à traiter 
aussi une question intéressante du maximum et du minimum du 
rayonnement sous une même latitude 4. Pour l’équateur, comme 
pour toutes les latitudes dans la zone torride, il y a évidemment 
2 maximum et 2 minimum du rayonnement qui sont donnés 
pour l’équateur dans le texte pol. Considérons maintenant la répar- 


et de Q donnés pour p=0,8 et p—0,7 dans le chapitre 


2.1047 4 


a dO à T 
tition de Qj ou -— pour les même © et pour les différentes 


latitudes; il faut remarquer que le maximum du rayonnement maura 
lieu que pour les latitudes qui ont le soleil au zćnite maisil est. de- 
placé vers la direction de jours croissants; on trouvera graphi- 
quement la position suivante du maximum: 
| ò= 80; Em = 12°30 

= 16°; Ym = 26030 
| Ô—=e ; ge =48°30 


(5) 


VIII. Quantité totale de la chaleur reçue par la terre. 


Dans le cas considéré le mouvement diurne et annuel 
de la terre nous n’intéresse pas puisque ces deux mouvements 
n’ont aucune influence sur la partie éclairée de la terre. Si 
W sera la quantité totale reçue par la terre on a: 


/ 1 
E ES | snn. p™.an 


Pour calculer l'intégrale de surface (1), nous utilisons les 
coordonnées géographiques A, +; on a alors: 


do = R? cos ọ . dẹ. dXA.sinh=cosę.cosh 


, dW _4CR2 f— 
Se 


f= sec$.sech cosżę, cosh . dà. dę 
o o 
de là: W =W' (0O—0% 
c-à—d. les quantités de chaleur sont proportionnelles aux accrois- 
sements de la longitude ©. 
Pour la quantité annuelle on aura: 


(3) were = 2 KR seo.psech Cor cosà. dÄ de 
0 0 
en divisant par Ach? on aura la quantité moyenne reçue par cm?. 
Pour p=l, c.—à—d à la limite de l'atmosphère on a: 


W = e E k cos . cosh. dh dę = 2R © 
o o 

et la quantité moyenne sera 

(4) | (P=1) QASET E= 269,5 kal. kg. 


Pour calculer Qu dans le cas général, considérons la fonction: 


1 
cos P.cosh 


(5) F (ę) = cosà. p „dh; 
J 
en posant: œ= secà . loge s 
; a 
et utilisant le développement comme pour e V1-¥ on arrive- 
ra au développement: 
F (g) = Yo + 1. sin?ę F ra . sinte +... 
Les coefficients y sont donnés par la formule (11) du texte pol. 
On aura ainsi pour W, annuelle et Qu annuelle la série suivante: 


AIR Ge, 15 S 
OZ ROLE LIGNE nt... 
> L E lus R Ca VE 57%, 4 | 
(4) d'Ae c | 9 lo ' an j 16 la SEITEN | 
On a ainsi: 


| p=0,8; Qu= 1,0716 ie = 183 cal. kg. 
(8) £ 


r 
| p=0,1; Qu= 0,8662 .— =148,6 cal. kg. 

puisque en veritć p est compris entre 0,8 et 0,7, on pourra dire 
que 38% de la chaleur est perdue dans l’atmosphère. 

IX. Applications de recherches spectrobolometriques. 


Comme nous montrent les études du spectre solaire, admis- 
sion du p constant pour tout le spectre, comme on a fait jus- 


qu'ici est inadmissible. Les pertes du rayonnement dans Pat- 
mosphère sont dues principalement à ła diffusion des rayons 
par les particules et à l'absorption de la chaleur par la vapeur 
d’eau. L'intensitć de la diffusion dépend de la longueur d'onde; 
l'air diffuse plus intensivement les vibrations bleu-violettes que 


les rouges et infra-rouges. Rayleigh a aussi donné une formule 


— H 


pour le coefficient de diffusion: K Sg où A est la longueur d'onde. 

La constante solaire c’est l'énergie totale du spectre à la limite 
de l’atmosphère; si donc ex. dh sera l'énergie des vibrations com- 
prises entre À et À+ dA, alors laire de la courbe représentative 
de eì sera égale à la constante solaire. De même si nous allons 
tracer la courbe ei représentative de l'énergie du spectre obser- 
vé à la terre, alors laire correspondante représentera l'énergie 
reçue par l’unité de surface. Les valeurs de ex nous allons pren- 
dre d'après les mesures de Abbot et Fowle, de même que les 
coefficients de transmission p pour les différentes longueurs d’onde 
(les valeurs suivantes de p sont valables seulement pour l'air sec): 


e RS E EE — 
À [00,30 |0p,35 01,40 Ou,45 0,50 | 0u,60 |0u,80| 14:,0 | Ip | 1,6 | 21,0 
10,539 | 2,684. 1,338 |6,027 6,062 5,042 | 2,665 1,657 | 0,898 | 0,532 | 0,247 
0,4 | 0,4 | 0,54 | 0,64 0,7 | 0,76 | 0,86 | 0,9 | 0,97 | 0,97 | 0,9 
| | I | 


A 
DA 


d’après ces valeurs on construira les courbes ex et éi pour le 
chemin traversć m; on a: 

(1) Eis DAD. E) 

Si on désigne par qo la constante solaire, par q l’intensité du rayon- 
nement, nous aurons: 


A 


| pa” .ea.da 


(2) Oz Go. À 


TEE 


Les figures correspondantes on trouvera dans le texte polonais. 

D'après la formule (2) on a trouvé une série des valeurs des 
composantes normales d'intensité du rayonnement qn (qa=q -sin h) 
pour les différentes valeurs de h (formule 9 du texte pol.) On a aussi 
développé q en série de Fourier suivante: 


(3) a On = 0,790 . sin h — 0,035. sin 3 h — 0,008. sin 5 h + 
0 


— 0,001 .sin 7 h—.... 

Mais dans (3) on a tenu compte seulement de la diffusion 
(l'air sec); le terme correspondant à l'absorption de la vapeur 
d'eau peut être pris empiriquement (d'après Angstró m) égal à: 

qo . 0,1.(0,85) M5 (fm) 0:275 
f— tension de la vapeur d’eau, à — densité de la diffusion. On 
aura le développement suivant: 
Sin h m %*%, 0,9" = 0,923 . Sin h 40,027. Sin 3 h + 0,005. sin 5 h + 
+ 0,001 KEE 
ce développement, joint avec (3), nous donnera l'expression sui- 
vante de la composante normale du rayonnement solaire, y com- 
pris l'absorption de la vapeur d’eau et de la diffusion: 


(4) e Qi = (0,790 — 0,092 . 0%) , Sin h — (0,035 +- 0,0027 
0 


„ [0275 ) X Sin 3 h + 
— (0,008 + 0,0005.1%*% ) Sin 5 h — (0,001 + 0,0001 . {%25 ) . sin7h—... 
Les valeurs calculées d'après (4) sont assez bien concordantes 
avec celles données expérimentalement. 
En se servant de la série de o on calculera les sommes 
de la chaleur pour les jours particuliers et de là pour les mois. 
Voici les résultats pour les valeurs mensuelles: 


Janvier 2,0 cal. kg. Juillet 16,4 cal. kg. 
Février 3,9 A Août 14,0 À 
Mars LOS Septembre 9,8 5 
Avril 12,0 $: Octobre 5,8 3 
Mai 15/65. 4, Novembre 2,9 

Juin 17,1 A Decembre 1,4 


„Prace* Wydz. III. M 15, Badania teoretyezne ilości ciepła. 


— 


An OS E 


pour la quantité annuelle Qr et pour Q;:et Q_e on aura: 
Q: = 108,4 cal. kg. 
Q+: = 85,6 cal. kg; Q—e= 22,8 cal. kg. 


Nous pouvons maintenant critiquer la prócision de la loi 
1 
= 


de Lambert; on doit en effet avoir p E '; mais appli- 


quant les valeurs précédentes de g on trouve p variable, notam- 
ment croissant avec m; la table suivante donne les valeurs de 


ki | que l’on pourra appeller les „coefficients moyens de trans- 
0 


mission“. 
h f=0m  f=2m  fz= gmm 
90° . . 0,818 0,710 0,656 
60° . . 0811 0,711 0,657 
40 . . 0,839 0,754 0,709 
30% . . 0,845 0,773 0,741 
20° . . 0,857 0,796 0,762 
10% . . 0,860 0,811 0,790 


La colonne f = 2%m est ordinairement appliquable en hiver et 
la colonne f=8"%m en été, puisque alors il y a plus de la va- 
peur d’eau dans l'atmosphère. Maintenant on peut bien expli- 
quer quelques inconvéniants qui résultent de l'application de la loi 


Q: 


de Lambert; par- exemple on trouve que le rapport Q est 


EE 
toujours trop fort pour les différentes valeurs de p; cela résulte de 
ce qwon a pris le même p en hiver qu’en été, malgré qu’en été 
le soleil est plus haut et est plus grand tandis qu’en hiver il en 
est au contraire et par conséquent, d‘après la table précédente, 
les valeurs de p en été sont inférieures à celles en hiver. On peut 
dire qu'en moyen il est plus naturel d'appliquer en été pour p 
les valeurs voisines de 0,7 et en hiver de 0,8. 
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